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1. QUARKS AND GLUONS 
1.1 Introduction 



In Standard Model the properties of quarks and leptons are remarkably similar, as far as the electroweak 
interactions are concerned. The quarks with six different flavours are grouped into the three doublets 
(u, d),(c, s),(t, b), in one-to-one correspondence to the three lepton doublets (u e , e),^, (i), (z/ T , r), as 
l/" - ) \ shown in the schematic chart in Fig. ^ Quarks and leptons interact in a similar way with the elec- 

troweak vector bosons 7, and Z. Furthermore, it is anticipated that the universal Higgs mechanism 
^ ■ is responsible for the generation of the quark and lepton masses (for more details see tTll 

In addition, quarks interact "on their own", revealing their specific property, the colour charge 
{colour), a conserved quantum number which is absent for leptons. A quark of a given flavour has three 
different colour states with equal masses and electroweak charges. Colour-induced interactions between 
quarks are mediated by gluons, the massless and electroweakly-neutral spin- 1 particles, 

In these lectures I will discuss the quark-gluon "corner" of the Standard Model, introducing the 
basics of Quantum Chromodynamics (QCD), the theory of quark-gluon interactions. Throughout this 
survey, the main emphasis will be put on the relation between QCD and hadrons, the observed bound 
states of quarks. This first lecture is devoted to the basic properties of the quark-gluon interactions. I will 
^ ■ frequently refer to Quantum Electrodynamics (QED), the more familiar theory of electromagnetic (e.m.) 

interactions, which is a useful prototype of QCD. 

In Fig. |3 the Feynman diagram of the electron-muon e.m. scattering is drawn together with the 
analogous diagram of the quark-quark interaction. For definiteness, I specify the quarks as having d and 
s flavours, the counterparts of e and \i in the Standard Model. The two interactions have many important 
similarities: 

- the colour-charged quarks emit and absorb gluons in the same way as the electrically charged 
leptons emit and absorb photons; 

- gluon and photon are massless; 

- both gluon and photon have spin 1 . 

Being guided by this analogy, one would expect that gluon exchanges generate a Coulomb-type 
interquark force, similar to the usual attraction/repulsion between the electrically charged particles. In 
reality, quark-gluon interactions are far more complicated. In particular, since the colour charge has 
three components, quarks can change their color states after emitting/absorbing gluons, as indicated on 
the diagram in Fig. 12 Hence, due to the colour conservation, gluons also carry colour quantum numbers, 
and, as a result, interact with each other. In fact, it is the gluon self-interaction that makes QCD-dynamics 
so peculiar. Yet QCD has one basic property in common with QED. Both theories have specific "internal" 
symmetries, named local gauge symmetries, to be discussed below. 




Fig. 2: Diagrams of lepton-photon (left) and quark-gluon (right) interactions. The dashed line is used to denote gluons and 
wavy line photon. 



1.2 Local gauge symmetry in QED and QCD 

The lepton-photon dynamics is described by one compact formula of the QED Lagrangian, 

L QED {x) = -- A F^{x)F^(x)+ Mx^iD^-m^ix), (1) 

l=e,fi,r 

where the physical degrees of freedom of the photon field (the electric and magnetic fields) are combined 
in the field strength tensor F^ u = d^A u — d u A^, with d v A^ = dj g^ x ^ ■ In the above, ipi(x) is the 
lepton Dirac field with spin 1/2 and mass mi, and a compact notation for the covariant derivative = 
+ ieA^ is introduced. Stalling from Lqed it is possible to derive Dirac equations for the leptons and 
Maxwell equations for the photon. Furthermore, Eq. Q yields the basics elements of the QED Feynman 
diagrams: the photon and lepton propagators and the lepton-photon interaction vertex. Employing these 
elements, one obtains the amplitude of a given e.m. process in a form of a perturbative expansion in the 
numerically small coupling a em = e 2 /47r. 

For us the most interesting property of Lqed is the gauge symmetry which reveals itself if one 
locally changes the phase of the lepton fields: 

ipi(x) ipKx) = exp [-1X0*0] iPl( x ) > 

Tpi(x) -> Tpi(x) =ijj l (x)exp[ix{x)] , (2) 

x(x) being an arbitrary function of the 4-coordinate x. If one simultaneously adds the derivative of the 
same function to the photon field: 

A„ -» Afa) = A„(x) + -d^xix) , (3) 

the Lagrangian © remains invariant. It is easy to check that the transformations © and (|3} leave intact 
the physical observables, such as the e.m. current "tpi^y^ipi or F^. Importantly, the local gauge symmetry 
implies that the photon mass vanishes, m 7 = 0. Indeed, adding to Lqed a mass term rn^A^A^ for 
the photon field would evidently violate the symmetry because the latter term changes under (|3j. The 
particular case \ = const (when only the lepton fields are transformed) corresponds to the global gauge 
symmetry which is responsible for the e.m. current conservation in QED. 

From the mathematical point of view the QED gauge transformations form a group. Let me re- 
mind that a given set of elements {gi} can be qualified as a group if three conditions are simultaneously 
fulfilled: 1) a multiplication rule can be defined g% * gk = gi, that is, a correspondence is established 
between a given pair of elements gi and g^ and a certain third element gi belonging to the same set; 2) 
the unit element go exists, so that go* gi = gi for each g^ and 3) the inverse element g^ 1 can be specified 
for each g^, so that gk * g^ 1 = go- All three conditions are valid for the (infinite and continuous) set 
of transformations © and © generated by the set of the arbitrary functions x( x )- Indeed, performing 
two gauge transformations one after the other, with Xi( x ) an d X2{x), is equivalent to the gauge trans- 
formation with Xvz{ x ) = Xi( x ) + X2{x). The unit element of this "multiplication " rule is simply the 
transformation with xo(^) = and the inverse element for each x(^) is — xi x ) ■ Importantly, the group 
multiplication is commutative (the group is Abelian) because the result of the overlap of two phase trans- 
formations is independent of their order. The group we are discussing is called U(l). Mathematically, 
it is equivalent to the group of rotations of the Cartesian coordinate system around one of its axes. The 
rotation angle plays the same role as the phase x- 

Gauge transformations in QCD have a more rich geometry. They are somewhat similar to the 
general rotations in the 3-dimensional space, involving 3x3 matrices, which do not commute. More 
specifically, a colour gauge transformation of the quark field 
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with a given flavour q = u,d, s, involves transitions between different color components, a sort of 
"rotations of colour coordinates": 



^ q (x)^^(x) = Ul(x)^ k q (x) 



(4) 



The elements of the 3x3 matrix Uj.(x) depend arbitrarily on the 4-point x. Furthermore, this matrix is 
unitary: 

U k = 5 k 3 , (5) 

or, in a symbolic form, U'U = 1. To explain why unitarity is necessary, let me invoke the following 
physical argument. Since photon does not distinguish quark colours, the only admittable form for the 
e.m. interaction of quarks is 



L em (x) = eQ q ^qk{x)l^q{x)A^(x) , 

fc=l,2,3 



(6) 



where Q q = +2/3 (—1/3) for q = u,c,t (d, s, b) and the summation goes over the quark colour indices. 
The gauge transformation © applied to the quark fields yields 



L em {x) - L> em (x) = ^ qi {x)U^{x^U{x))^{x)A^{x) . 

Clearly, only if © is valid, L em remains invariant. 

The usual exponential representation of the gauge transformation matrix is: 



(V) 



U l k (x) = exp 



(A a )l 



a=l 



(8) 



It contains eight independent and arbitrary functions x a ( x ) multiplied by eight reference matrices A a 
(a = 1, ...8). The latter have the form chosen by Gell-Mann: 
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(9) 



The A-matrices have the following properties: A at = A a (hermiticity), Tr\ a = 0, and [A a , A 6 ] = f abc \ c 
(noncommutativity), where f abc are totally antisymmetric constants (/ 123 = — f 213 , etc.). It is easy to 
check that U (x) defined in (|3l obeys unitarity and has a unit determinant det U = 1. The (infinite and 
continuous) set of noncommutative U matrices forms a group, which is called SU (3). One may ask: why 
there are eight independent functions x a determining the rotations of the three color states? Isn't eight 
too many? The point is that the quark fields ip l are complex functions, hence the matrix U(x) is also a 
complex function of x determined by 2 x 9 real functions. Only 8 of them are independent because there 
are altogether 10 constraints: nine provided by the unitarity relation © and one by the unit determinant. 
With eight "rotation angles" the group SU (3) is of course quite different from the group of rotations in 
three dimensions which has only 3 parameters. 



1.3 QCD Lagrangian 

The QCD Lagrangian has to be exactly symmetric with respect to the local gauge transformations © 
with matrices ©. This property serves as a guiding principle for constructing Lqqd- 

We start from the part of the Lagrangian which describes the propagation of free quarks: 

L qU ark{x)= ^ 12 A ^(id^ - TUq)^ (x) \ , (10) 

q=u,d,s,... \fc=l,2,3 / 

and yields the usual Dirac equation for spin 1/2 particle for each quark with a given flavour and colour. 
Similar to the case of QED, the expression dTOb is not invariant with respect to the local gauge transfor- 
mations ©. An additional term with derivatives of x a {x) remains: 

Lquark{x) ^ L quark { X ) + £ % i (x) [iU\ * (x)d fl U} (x)] -flft (x) . (11) 

q=u,d,s,... 

To restore gauge invariance one follows the same scenario as in QED or in the electroweak theory 
The idea (put forward long ago by Yang and Mills) is to introduce "compensating" spin-1 fields interact- 
ing with quark fields. There should be one separate spin-1 field for each of the eight degrees of freedom 
determining the gauge transformations ©. In this way eight gluons enter the game, with the following 
quark-gluon interaction term added to the Lagrangian: 

L in t(x)=g s hiW^-ftftoWx). (12) 

q=u,d,s,... 

Here g s is the dimensionless coupling analogous to e in Lqed- In contrast to e.m. interaction, where 
the photon field is electrically neutral, the gluon fields also carry colour charge, so that the colour state 
of a quark changes after emitting/absorbing a gluon. The colour of the gluon distinguished by the index 
a = 1...8 can be identified with a superposition of quark and antiquark colours. For example the gluon 
field A l Ax) is in the same colour state as the quark-antiquark pair ^1^ + ^52^- The A^-fields in (fT2l> 
have to be gauge-transformed in the following way: 

^A«(x) - U(x)^A;(x)U\x) - -d,U(x)U\x) , (13) 
z z g s 

so that the overall change of Li nt cancels the symmetry breaking term on r.h.s. of (fTTT l. It is a simple 
exercise to check that the combination of transformations © and (fT3l leaves the sum Lg Uar A; + Li nt in- 
variant. We see that in QCD the "compensating" transformation of gluon fields (fT3l is more complicated 
than its analog © in QED: the addition of the derivative over x-functions is accompanied by a "color 
rotation". 

To complete the Lagrangian one has to add a gauge-invariant term describing the propagation of 
gluon fields : 

L glue (x) = -±G%(x)G a ^(x), (14) 

where 

G% = d,Al - d v A% + g s P hc A\Al , (15) 

is the gluon field-strength tensor. The local gauge invariance of L g i ue implies that gluons are massless. 
At the same time, G® u is a more complicated object than its QED analog F^ v . Indeed, substituting 
IT31 to (fT4l . we notice that not only the terms quadratic in A" emerge (propagation of gluons) but 
also the three- and four-gluon vertices (gluon self-interactions). Note that formally, both properties of 



Quark propagator 


<0|T{^(s)^ fc (0)}|0> = iS{ J d'p e"^ 


Gluon propagator 


(0\T{Al(x)A b u (0)}\0) = -i5 ab J d 4 k^-e ikx 


Quark-gluon vertex 




3-gluon vertex 


-f f abc [d^A a u (x) - d v Al{x)]A b ^{x)A cv {x) 


4-gluon vertex 


- 9 ^f abc fadeA b ll (x)Al(x)A^(x)A^(x) 



Table 1: Propagators and vertices in QCD. 



gluons: colour quantum number and self-interactions are due to the noncommutativity of the gauge- 
transformation group (the fact that f a ^ c / 0). The final form of QCD Lagrangian is obtained by adding 
together the three pieces introduced above: 

LqCD = L g i ue + Lq U aTk + ^int 

= -\g%G^ + 4>SD»r - m q )^ q , (16) 

where = — ig s ^-A a ^. To summarize, Lqcd describes not only quark-gluon interactions but also 
gluodynamics, the specific gluon self-interactions which have no analog in QED 1 . 

The quark and gluon propagators and vertices derived from Lqcd are collected in Table 1. The 
formula of the gluon propagator has a certain degree of freedom related to the fact that the physical 
massless gluon has only two polarization/spin states whereas the field A® has four components. To make 
things working, one follows the same recipe as in QED. An additional constraint on the gluon field is 
introduced, the so called gauge-fixing condition. The gluon propagator given in Table ^corresponds to 
the usual Feynman gauge. Note that in QCD purely gluonic loop diagrams are possible, in which case 
one has to take care of subtracting the contributions of unphysical components of A" also in these loops. 
It is usually done by adding specially designed fictitious particles (the so called Fadeev-Popov ghosts) 
which only appear in the loops, and are not shown in Table 1. 

Feynman diagrams in QCD are obtained employing the quark-gluon propagators and vertices as 
building blocks. However, the use of diagrams makes sense only if the perturbative expansion in g s is 
meaningful. To obey this condition, the quark-gluon coupling 

« s = fS (17) 

4-7T 

'in QCD "light emits light" at the level of the fundamental interactions entering Lagrangian, as opposed to QED where light- 
by-light interaction appears only as 0{ptl m ) quantum correction (when photons exchange virtual electrons via box diagrams). 
For a classical Maxwellian electrodynamics, self-interacting e.m. fields would mean "new physics beyond standard theory". I 
am not aware of any discussions of photon self-interactions in the times before quantum field theory. Interestingly, the light 
emitting light was mentioned in poetry. I found the following sentence written in XIII century by Armenian poet Kostandin 
Erznkazi |2|: "And so the light was born from the light, the great light of Sun..." (in translation from Armenian) 



Fig. 3: Some of the diagrams corresponding to quantum loop corrections to the quark-quark scattering 



the QCD analog of the e.m. coupling a em = e 2 /At:, has to be sufficiently small, a s <C 1. If this 
condition is fulfilled, then, for example, the 0(a s ) diagram of the quark-quark interaction in Fig. |2] 
dominates over the higher-order diagrams with additional gluon exchanges. One has then a tractable 
situation, with quarks and gluons propagating quasi-freely. Thus, there is an important question to be 
addressed: how large is a s ? 

1.4 Running of the quark-gluon coupling 

To answer the above question, one has to investigate quantum effects in QCD, i.e., creation and annihi- 
lation of virtual gluons and quarks described by loop diagrams. In QED, which we use as a prototype, 
quantum loops generate very small effects, such as Lamb shift (the correction to the Coulomb force due 
to the virtual electron-positron pairs) or the 0(a em ) one -loop correction to the muon magnetic moment. 
These effects, being accessible in precision experiments, play a minor role in the bulk of electromagnetic 
processes. 

In QCD, quark-gluon loops are far more pronounced and play a crucial role in determining a s . To 
have a closer look, let us consider the quark-quark scattering amplitude. In the lowest-order (at the tree 
level) the amplitude is given by the one-gluon exchange diagram shown in Fig.|2]for the particular choice 
of quark flavours. Having at hand the Feynman rules in QCD it is easy to write down the amplitude: 



where q 2 < is the momentum transfer squared. Below I will also use the notation for the momentum 
scale: Q = \J~Cp, where Q 2 = —q 2 > 0. Considering 0(a s ) corrections to the amplitude ( fT8t one 
encounters diagrams shown in Fig. [3] They contain gluon or quark loops inserted within the exchanged 
gluon line or in the vertices. These loop effects turn out to be extremely important for evaluating a s . 

Let me first explain the loop diagram calculation in more detail. After substituting propagators 
of virtual particles, one arrives at four-dimensional Feynman integrals over the 4-momentum k flowing 




(18) 



inside the loop, typically: 



Iloop{q,m q 



d 4 k 



(2vr) 4 (k 2 - m 2 )({k + q) 2 - m 2 ) ' 



(19) 



where the explicit expression in the numerator (depending on the spin of the loop particles) does not play 
role in our discussion. For simplicity, I will also put to zero the quark masses, a reasonable approximation 
if Q S> m q . To calculate loop integrals, one usually employs the method of dimensional regularization. 
The idea is to lower the number of dimensions in the integral (fl9l making it convergent. One replaces 
4 by a generic integer number D, then calculates the integral as a function of D and after that considers 
the result at an arbitrary noninteger D = 4 — e, schematically: 



Iloo P (q, 0) = J d 4 kf(k, q) - // 4 - D ) J d D kf(k, q) = I(q, (i, D) - I(q, M , 4 - e) 



(20) 



The auxiliary mass scale fi is introduced to keep unchanged the physical dimension of the integral. The 
major advantage of dimensional regularization is in preserving gauge symmetry of the amplitudes at each 
step of the calculation. The divergent part of the integral in this scheme is represented in a form of terms 
proportional to 1/e. A generic expression for the loop integral has a form: 



iioo P (q,o) 



A-D 



dPk 



(2tt) d k 2 (k + qf 



h \og{-q 2 /^) + -+I 2 
e 



(21) 



where are calculable finite coefficients. 

QCD is a renormalizable theory (similar to QED), which means one can absorb all divergent 
1/e terms into the so called Z -factors. Multiplying by Z~ l the coupling g s , and by the corresponding 
factors the quark masses, quark and gluon fields, one defines the finite (renormalized) quantities; e.g., the 
renormalized coupling is g r s en = Zg X g s . Having in mind the validity of the renormalization procedure I 
will simply ignore divergent terms appearing in d2TT i and in other loop integrals. 

Adding the 0{a s ) diagrams in Fig. [3] to the tree-level amplitude, results in the same expression 
(1181 . but with a s replaced by an effective coupling depending on the momentum scale: 
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where a shorthand notation 
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(22) 



(23) 



is introduced. Here rif is the number of "active" quark flavours in the loop diagrams. Only those quarks 
which have m q <C Q contribute to a e /f (Q 2 ). Importantly, 0o is positive, because the term 1 1 originating 
from the gluon loops exceeds the quark-loop contribution 2nj/3 (since nj < 6 in any case). 

Taking ay* at a different scale Qo, 
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(24) 



and dividing d22l by (124b one obtains, with an accuracy of 0(a 

v e //(Q ) 



2): 
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(25) 



the relation between effective couplings at two different scales. The approximation d25l is valid only if 
ay* (Qo) is sufficiently small and the higher-order corrections are negligible. Remarkably, d25b predicts 



a s logQ + ( a s logQ ) + 

Fig. 4: One- and two-loop diagrams contributing to the running ofa s . 



that a e J ' (Q) decreases when the momentum scale Q increases. Thus, if the perturbative expansion in 
a s is applicable at certain Qq, it should behave even better at Q > Qq. 

Note that ( I25t still has to be improved. At Q — > oo the logarithm becomes very large driving 
the combination a s log Q 2 to rather big values, so that the O ([a s log (Q 2 /Qo] 2 ^j correction originating 
from the two-loop diagrams shown in Fig.|4]becomes important, and the whole perturbative construction 
is again in danger. Fortunately, a systematic resummation of all O ({a s log (Q 2 /Q 2 ,] n ^j corrections is 
possible. In practice, one does not need to calculate all multiloop diagrams, which would be a tremendous 
work. Instead, the renormalization-group method is used, which is however beyond our scope. The 
resulting expression for the running coupling is: 

asiQ) = ~ MQo)o Q 2 ■ (26) 



l+^22/3 .lo, 



-0 



(hereafter the superscript 'eff is omitted). Expanding the denominator in d^ol i and retaining only the 
first two terms we return to the relation d25t . 



1.5 Asymptotic freedom 

A consistent use of the running QCD coupling is possible if one can find a range of Q where a s (Q) 
is numerically small. The first indications that a s is indeed small at large momentum transfers were 
obtained in the studies of deep-inelastic lepton-nucleon scattering (to be discussed in Lecture 3). This 
remarkable discovery paved the way for using QCD perturbation theory with the running coupling in 
many other processes. 

To trace the numerical behavior of a s (Q) one has to fix the coupling at a certain large scale us- 
ing an experimental input. One possibility is the decay of Z-boson to a quark-antiquark pair. Quarks 
originating in this decay inevitably build hadrons in the final state (see the next subsection). To avoid 
hadronic complexity, one measures the total inclusive width T(Z — ► hadrons), so that the probabilities 
of all possible quarks — > hadrons transitions add up to a unit. The majority of Z — > hadrons events ob- 
served at LEP has a spectacular structure of two hadronic jets originating from the initial, very energetic 
quark and antiquark (E q = Eq = mz/2 in the Z rest-frame), see, e.g. [3 1. On the other hand, the share 
of > 3-jet events in Z — > hadrons, with additional jets originating from gluons and/or from secondary 
quark-antiquark pairs, is small. This observation clearly indicates that the quark-gluon coupling at the 
scale mz is small, or in other words the initial quark pair interacts weakly during the short time after its 
creation. 

The perturbative diagrams of Z — > qq ( q = u,d, s, c, b) including the gluon emission Z — > qqG 
are shown in Fig. |5] Evaluating these diagrams one gets the perturbative expansion for the total hadronic 
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Fig. 5: The lowest-order diagram (a) and some of the higher-order in a a diagrams (b-d) determining the total width Z 
hadrons. 



width, schematically: 

T(Z — > hadrons) = 



q=u,d,s,.. L 



T(Z -> qq) (l + Cla s (m Z ) + C q 2 a 2 s {m z 
+T(Z - qqG) (l + Cfa s (m z ) + ...) + ... 



(27) 



where T(Z — > gg) and T(Z — > gf/G) (the latter starting from 0(a s )) are the corresponding perturbative 
widths, and Cf 2 ,Cf G ,.. are the calculable coefficients. Smallness of a s allows to neglect all higher- 
order corrections starting, say from 0(a~). In the above expression a s is taken at the scale mz, the 
characteristic scale in this process. One can trace how the running coupling builds up in (127b . by summing 
up all logarithmic corrections generated by the loop insertions similar to the one shown in Fig. |4j;. 

Comparing the result d27b of the theoretical calculation with the experimental data yields [ 4 1 : 



a s (m z ) ~ 0.12. 



(28) 



As expected, it is a rather small number, so that the whole perturbative treatment turns out to be quite 
consistent. Using d26l and d28t . one predicts a s (Q) at lower scales. The curve plotted in Fig. [6]is taken 
from [4] and reflects the current status of the running coupling, including all known (and reasonably 
small) higher-order corrections to d26b . As can be seen in this figure, there is a wide region spreading up 
to Q ~ 1 GeV, where a s is small and perturbative QCD is applicable. Furthermore, a s extracted from 
various processes at different scales agrees with the running behavior predicted in QCD. The most spec- 
tacular consequence of d26b is the vanishing of the running quark-gluon coupling at Q — > do, revealing 
that QCD is asymptotically free. 



1.6 Confinement and hadrons 

Quite an opposite situation takes place in the quark-gluon interactions at small momentum transfers (at 
long distances). According to d26b . if one starts from Q » 1 GeV and drifts towards smaller scales, a s 
grows to O(l) at Q < 1 GeV (see Fig. 0. Perturbation theory becomes useless in this region, because 
e.g. in the quark-quark scattering, an infinite amount of higher-order quark-gluon diagrams has to be 
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Fig. 6: The dependence ofa s (jJ,) on the running scale p [41: the upper and lower curves indicate the theoretical uncertainty, 
the points with errors are the values of a s extracted at different scales using various methods. 



taken into account. Moreover, at a certain momentum scale denoted by Aqcd the denominator in d26t 
vanishes and a s (Q) diverges. The relation between a s {Q) and Aqcd following from d26l is: 



I skip some details concerning the dependence of Aq^ d on the number of active quark flavours, and on 
the theoretical scheme of the loop-diagram evaluation, relevant for the higher-order corrections to ( 129 i . 
The important fact is that d29b and d30t were derived neglecting the quark masses. We could have put 
to zero the quark masses in Lqcd from the very beginning, starting with a theory without dimensionful 
parameters. The emergence of an intrinsic energy scale in a theory with dimensionless coupling g s 
{dimensional transmutation) is a specific property of QCD, due entirely to quantum effects. 

The breakdown of perturbation theory and the exploding behavior of a s (Q) at Q — > Aqcd are 
actually anticipated. Long before QCD was invented it was known that at long distances quarks and 
antiquarks strongly interact and form hadrons, the quark-antiquark (meson) and 3-quark (baryori) bound 
states. The properties of hadrons will be considered in a more systematic way in the next Lecture. For 
the present discussion it is important that the characteristic energy scale of hadronic interactions is of 
O(Aqcd)- Hence, it is quite natural that the formation of hadrons is due to the strong, nonperturbative 
quark-gluon force emerging in QCD at momenta ~ Aqcd- 

Moreover, hadronic matter is the only observable form of quarks and gluons at long distances. 
In any process of quark and antiquark production, independent of the energy/momenta involved, quarks 
form hadrons in the final state 2 . Note that in QED the situation is quite different: isolated leptons 
and other electrically charged particles are observed, and the e.m. bound states (e.g., hydrogen atom, 
positronium or muonium) can always be split into constituents if a sufficient energy is supplied. 

2 There is one exception: t-quark, decaying via weak interactions, is too short-lived to be bound by quark-gluon forces. 
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The experimentally measured value d28b corresponds, roughly, to 



Aqcd = 200 - 300 MeV. 



(30) 
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Fig. 7: Schematic view of the a s behavior at different scales. 



In QCD, the non-observation of free colour-charged particles (quarks, antiquarks and gluons) is 
arranged in a form of the colour confinement principle, postulating that all observable states, i.e., all 
hadrons, have to be colour-neutral. The confinement principle was never rigorously proved, because of 
our limited ability to work with QCD beyond perturbation theory. Nevertheless, all experimental results 
concerning hadrons, as well as lattice simulations of QCD at long distances, unambiguously support 
colour confinement. 

Before the era of QCD the search for free quarks was very popular among experimentalists. The 
fractional electric charge was a smoking-gun signal to be observed. The hunt for quarks was a part of 
many accelerator experiments. Not surprisingly, free quarks have never been found at accelerators or in 
other places (in cosmic rays, water, ice, meteorites etc.). Nowadays, one would hardly invest efforts in 
the search for free quarks. We are confident that QCD obeys confinement. 

1.7 Quark masses 

The current intervals of quark masses presented in the Particle Data Tables [4| are: 

m u = 1.5 -T- 4.5 MeV, m c = 1.0 -5-1.4 GeV, m t = 174.3 ± 5.1 GeV; 
m d = 5 -T- 8.5 MeV, m s = 80 + 155 MeV, m b = 4.0 + 4.5 GeV . (31) 

The fact that the masses are spread over five orders of magnitude, is a reflection of some fundamental 
flavour physics not related to QCD, e.g. the Higgs mechanism of the Standard Model. Hence, the 
masses m q entering QCD Lagrangian (fT6t are "external" parameters. At the same time quark masses 
are evidently changed in the presence of quark-gluon interactions. At long distances, within hadrons, 
each quark acquires, roughly speaking, an extra addition of O(Aqcd) to its "bare" mass. It is however 
very difficult, if not impossible at all, to define this constituent quark mass in a model-independent way. 
Representing the hadron mass as a sum of the constituent quark masses plus some interaction energy, 
e.g. for a meson: 

M meson = m c ° nstit + mf nstit + E int , (32) 

one can always redistribute the part of Ei nt between the quark and antiquark masses. Since free, on-shell 
quarks are not observed, the usual definition of the particle mass (the minimal possible energy of the 
one-particle state) does not work. 

The mass values presented in JTil have nevertheless quite a definite meaning. These are "short- 
distance" masses of the virtual quarks. As we already know, at large virtualities quark propagation is 



(a) 



(b) 



Fig. 8: The free quark propagator (a) and the gluon correction to it (b). 



quasi-free and a consistent use of perturbation theory derived from QCD Lagranigian is possible. In 
particular, the free quark propagator 

S(p) = ^±^i (33) 
p A — m q 

is applicable at \p 2 \ ^> Aq CD , m 2 , with the bare mass m q from Lqcd- Furthermore, at short distances 
the quark-gluon interactions are calculable in terms of series in a s . The quark propagating at short 
distances can emit and absorb a gluon (see diagram (b) in Fig . |8j . From the experience with the running 
coupling one expects such loops to be important. Adding the loop diagram to the free propagator yields 
the same expression as d33l with m q replaced by an effective mass depending on the momentum scale 



K ff (Q) 



1 - 7— ( 70 log ^2 + const 



(34) 



where I omit the divergent part, knowing that it will be absorbed by renormalization. In this case the scale 
Q ~ VTp^I i s determined by the virtuality of the quark and 70 = 4 is the result of the explicit calculation. 
Again, as in the case of a s , we can relate the effective masses at two different scales. Writing down the 
above equation for another scale Qo an d using instead of m e q ^ a conventional notation m q we obtain: 

m q (Q) = m q (Q ) (1 - (I) ^M/? log |H , (35) 

where I have multiplied and divided the logarithmic term by (3q and used a s = a s (Qo) which is correct 
with 0(a s ) accuracy. With the same accuracy the expression in parentheses can be transformed further: 

7oWQo) fll Q 2 (, a s (Q ) fll (f S 
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rn q (Q) = m q (Q ) ( , (37) 



Using d25l we obtain 

rh.JO\ = rh.JO^ ( 

the formula for the running mass. A more rigorous derivation is possible using the renormalization group 
method. Also in recent years, the higher-order corrections to d37b have been calculated. 

The masses presented in [4 1 are the running masses (also called MS masses if one specifies the 
particular renormalization procedure) normalized at some large scale. The light u,d,s quark masses are 
traditionally taken at Q = 2 GeV, whereas a more appropriate scale for the heavy c and b quark masses 
is the quark mass itself, Q = m c and Q = mi,, respectively, (which means, e.g., the virtuality of the c 
quark is p 2 = —m 2 ). The fact that quark masses run with the scale, is in accordance with the absence of 
isolated quarks among observable states. 



1.8 Two branches of QCD 

To summarize, QCD yields two qualitatively different pictures of quark-gluon interactions: 



1) at high momentum-transfers i.e., at short average distances, perturbative expansions in a s are 
applicable in terms of Feynman diagrams with quark and gluon propagators and vertices. In this region 
the scale-dependence (running) of the coupling and quark masses should be properly taken into account. 

2) at low scales, that is, at long distances, one loses control over perturbative interactions between 
individual quarks and gluons; the latter strongly interact and form hadrons. 

Accordingly, QCD is being developed in two different directions. The first one deals with short- 
distance physics accessible at high energies. One studies specific processes/observables calculable (at 
least partly) in a form of a perturbative expansion in a s . A typical short-distance process is the jet 
production in Z decays considered above, other examples will be presented in Lecture 3. 

The second direction deals with nonperturbative quark-gluon interactions at long distances and 
with hadron dynamics. A complete analytical evaluation of hadronic masses and other parameters di- 
rectly from Lqcd is not yet accessible. Instead, a powerful numerical method of simulating QCD on 
the space-time lattice has been developed. Lattice QCD has become a separate field, which is beyond 
the scope of these lectures (for a pedagogical introduction see e.g.,[5|). Still there is a lot of interesting 
advances in the long-distance 'branch" of QCD, so that I will only be able to cover a part of them. As 
demonstrated in Lecture 2, many important features of hadron spectroscopy follow from QCD at the 
qualitative level. The relation of long-distance dynamics to the nontrivial structure of the QCD vacuum 
will be discussed in Lecture 4. An approximate analytical method of QCD sum rules based on this 
relation and used to calculate hadronic parameters will be overviewed in Lecture 5. 

One might think that physics of hadrons plays a secondary role, because the most important di- 
rect tests of QCD in terms of quarks and gluons are done at short distances. Let me emphasize the 
fundamental importance of hadron dynamics by mentioning two topical problems: 

1) The origin of the nucleon mass 

Proton and neutron are the lowest and most stable baryons, with the quark content uud and udd, 
respectively. Their masses 

m p ~ m n ~ 940 MeV , (38) 

are substantially larger than the tripled quark mass m u ^ = 0{ few MeV). We conclude that ~ 99% 
of the baryon matter in the Universe is not related to the "fundamental" quark masses generated by the 
Higgs mechanism or from some other flavour dynamics. The bulk of the baryon mass is due to the 
long-distance quark-gluon interactions. Indeed, adding a "constituent" mass of O(Aqcd) to each quark, 
one gets the order-of-magnitude value of m P)n . Certainly, the problem of the nucleon mass is quite 
fundamental and has to be solved within long-distance QCD. 

2) Extracting electroweak parameters from B decays 

The weak decays of B mesons (the bound states of b quark and light antiquark) represent a valu- 
able source of information on fundamental aspects of electroweak interactions, such as the quark-mixing 
CKM matrix and the origin of CP-violation (see [6]). One topical example is the B — > -kIvi decay ob- 
served at B factories EJ. This decay (see Fig.|9j is driven by the b — > u weak transition, proportional to 
Vub, one of the poorly known CKM matrix elements. In order to extract this fundamental parameter from 
the experimentally measured partial width, one needs to divide out the hadronic B — > it transition am- 
plitude (form factor). The latter is essentially determined by the long-distance quark-gluon interactions. 

Below, in Lecture 5 we will discuss the (approximate) solutions of the two abovementioned prob- 
lems. 



Fig. 9: B — > tyIvi decay involves a hadronic transition between B and ty states. 



2. FROM QUARKS TO HADRONS 
2.1 Mesons and baryons 

Let us now have a closer look at the properties of hadrons 3 , the bound states of quarks. There are lots 
of experimental data on hadrons accumulated in the Review of Particle Physics [4|. As we shall see, 
employing various symmetries of QCD, it is possible to predict, at least qualitatively, many observable 
regularities of hadronic spectra and interactions. 

The pre-QCD quark model of hadrons failed to explain why only mesons (qq), baryons (qqq) and 
antibaryons (qqq) are observed. Why, for example, the diquark (qq) or four-quark (qqqq) bound states 
are absent? In QCD, one immediately gets an explanation based on the colour confinement principle. 
The quark-antiquark and three-quark combinations can form colour-singlet states, whereas diquark and 
four-quark compounds are always colour-charged. The colour-singlet meson state is obtained by sim- 
ply summing over the colour indices of the quark and antiquark. For example, the 7r + -meson has the 
following flavour/colour structure: 



The way baryons are built is less trivial. Two quarks are arranged in a coloured diquark state eijk\q^q k ), 
where is totally antisymmetric with ei23 = 1. This state has the colour SU(3)-transformation prop- 
erties of an antiquark (with a colour index i). Hence, we can obtain a colour-neutral state combining 
diquark with the third quark and summing over colours in the same way as in mesons. One important 
example of a three-quark baryon state 4 is Q~ baryon with spin 3/2: 



With respect to the flavour and spin, this state is totally symmetric: all quark have the same flavour and 
parallel spins. Therefore, colour degrees of freedom provide the antisymmetry demanded by the Pauli 
principle for any bound state of fermions. 

3 The word 'hadron' was coined by Okun in |7| where he wrote: "It is reasonable to call strongly interacting particles 
hadrons, and the corresponding decays - hadronic. In Greek the word hadros means "large" "massive", in contrast to the word 
"leptos", which means "small" "light". The term hadron refers to long-lived mesons and baryons, as well as to resonances." 

4 Note that if there were four colours of quarks, with the corresponding SU(4) symmetry, baryons would have been built 
from four quarks, with profound consequences for the physical world, e.g., atoms with fractional electric charge. 



(39) 




(40) 



2.2 Quark model of hadrons 

The hadron decompositions d39l and d40b resemble chemical formulae, displaying the content of a com- 
posite state. In reality hadrons are far more complicated than atoms and molecules, because the masses 
of the u, d, s quarks are smaller than the QCD long-distance scale: m Ui d -C Aqcd and m s < Aqcd- 
Hence, light quarks are purely relativistic and the number of quarks, antiquarks and gluons within a 
hadron cannot be fixed. Since QCD is a quantum field theory, additional quark-antiquark pairs or gluons 
are created and annihilated inside the bound states virtually, i.e. within the time/distance intervals of 
0(1/ Aqcd)- As a result, the general decomposition of the physical pion state is not simply d39l but 
rather 

I 7T+) = \u%) + I ui diq 3 qj)+ \u k {\ a )\d t G a ) + ... , (41) 

q=u,d,s,.. 

where G denotes a gluon, the Lorentz indices are not specified and a summation over colour indices 
is implied. In the above sum, the first term represents a state with the minimal particle content (the 
so called valence-quark content), and ellipses indicate all other multiparticle fluctuations. Naturally, all 
components of the pion state have to be colourless, with the same ud overall flavour, to obey the colour- 
neutrality and flavour conservation. We conclude that hadrons are, in general, many-body systems with 
relativistic constituents. Therefore, simple quantum-mechanical models with an interquark potential 
cannot fully describe pion or other light-quark hadrons. 

Before attempting to solve the QCD dynamics, it is useful to apply the symmetries of QCD La- 
grangian. The space-time (Lorentz-Poincare) invariance implies that the total angular momentum (or 
total spin) J of a hadron is a well defined and conserved quantum number. In addition P- and C- parities 
are conserved in QCD (as opposed to the electroweak theory). Therefore, for a given hadron, the spin- 
parity combination J p (J PC for flavour-neutral hadrons) is the next important signature after the mass. 
Spin-parities are indicated for each observed hadron in its entry in |4j. 

To proceed in hadron spectroscopy, let us have a more detailed look at the mesons having the same 
flavour content ud as 7r + . Each meson state is a complicated coherent decomposition similar to (14 li t. 
Nevertheless, since J p is conserved, it is sufficient to consider the valence-quark component to count 
all possible combinations of J p starting from the lowest possible spin. The total angular momentum 
of the valence quark-antiquark state is a sum of the quark and antiquark spins plus the orbital angular 
momentum: 

J = S + L, (42) 

where 5 = s q + Sq is the total quark spin; 5 = 1(0) if the individual spins are parallel (antiparallel). 
Accordingly, there are two possible states with L = and J = S. One is with J p = 0~ (pion) and 
the other with J p = 1~ (p meson). The negative P-parity attributed to these states is obtained from the 
following rule: P = — (— 1) L . Notice an additional minus which has to be added to account for the so 
called "intrinsic" P-parity of the relativistic quark-antiquark system. Turning to the states with L = 1, 
one encounters three mesons with 5 = 1 (J p = + , 1 + , 2 + ) and one with 5 = (J p = 1 + ). They are 
listed in Table |2] according to the classification of [4]. A similar counting can be done for L = 2, 3, 
predicting J > 2 mesons. Some of them can be found in |4). Generally, it is rather difficult to observe 
hadrons with higher spins. Having larger masses, these states have many decay channels and, therefore, 
a large total width, complicating their experimental identification in a form of a resonance. Baryons 
from u, d, s quarks with different J p listed in [4| can also be interpreted, at least qualitatively, in terms 
of three-quark valence states with the orbital momentum L between diquark and the third quark. 

The angular momentum is not the only source of excited hadron resonances. There are mesons 
which have the same J p as 7r or p but a larger mass. These states are somewhat similar to the radially 
excited levels in a potential. For the pion a natural candidate of such excitation is the 7r'(1300) state 
with J p = 0~, whereas p meson has at least two experimentally established radially-excited partners 
with J p = 1~: //(1450) and p"(1700) [4|. Note that because the orbital momentum L is not conserved 
in relativistic theory, the L = 2 state with J p = l - cannot be simply distinguished from the "radial 
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Table 2: Spectrum of the lowest ud states. 



excitation" of the L = state with the same spin-parity, - another difficulty for the potential models of 
light-quark hadrons. Ultimately, one has to think in terms of purely relativistic extended objects, some 
kind of quark-gluon strings having a spectrum of "radial" and J excitations. However, attempts to derive 
a string picture for hadrons directly from Lqcd were not successful so far. 

Quark-gluon interaction is flavour-independent. Therefore, given that a ud meson with a certain 
J p exists, the mesons with the same J p containing all possible quark-antiquark flavour combinations 
should also be observed. The flavour partners of the pion (p-meson) with J p = 0~ ( J p = 1~) are listed 
in Table|3] Almost all of them have been observed; the masses and other characteristics are given in |4j. 
The only temporary exceptions are the pseudoscalar bb meson (%) and the vector be meson (B*). These 
two states are not yet in [4|, due mainly to experimental reasons. 

The heavy quarkonia, i.e., the mesons consisting of heavy quark and antiquark (cc or bb) are of a 
special interest. Here the masses of interacting quarks are large enough compared to their characteristic 
energies within hadrons: 3> Aqcd- In other words, heavy quarks are nonrelativistic objects with 
respect to QCD interactions. It is therefore possible to approximate the quark-gluon interactions with a 
nonrelativistic potential, putting the hadronic calculus on the safe ground of quantum mechanics. The 
Coulomb quark-antiquark potential V(r) = a s /r, derived from the one-gluon exchange, is valid at small 
distances. To provide quark confinement, a certain long-distance part of the potential, infinitely growing 
at r — > oo should also exist (e.g., oscillator or linear potential). This part of the potential cannot be 
directly calculated from Lqcd and is usually modeled and fitted to the observed quarkonium spectra. 
Importantly, numerical studies of QCD on the lattice confirm the existence of the confining potential 
force between heavy quark and antiquark. 

2.3 Isospin 

In addition to the exact colour- and space-time symmetries, QCD Lagrangian possesses approximate 
flavour symmetries originating from the pattern of quark masses. Since the latter are generated by some 
external mechanism, the flavour symmetries do not have fundamental roots in QCD. Nevertheless, they 
provide very important relations for hadron masses and hadronic amplitudes. 

Let us start with the u and d quarks and rewrite the QCD Lagrangian, isolating these two flavours: 

Lqcd = Tpu(iD^ ~ m u )ip u + $ d (iD M y - m d )ip d + L glue + L SjC)6)t . (43) 
The smallness of the u and d masses, m u ~ <C Aqcd, implies that also their difference is small: 

m d -m u <. Aqcd ■ (44) 
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Table 3: Pseudoscalar ( J p =0 ) ( upper lines) and vector ( J p = ) (lower lines) mesons with different flavour content. 



Neglecting this difference we use a new notation for the common u, d quark mass: 

m u ~ rrid ~ rh . 

In this approximation, 

Lqcd * L^cd = ~*{D^ L -m)m + L glue + .. , 
where a new, two-component spinor field (doublet) is introduced: 



(45) 



(46) 



The theory described by the r.h.s. of d46t is not exactly QCD, but is very close to it. The new Lagrangian 
Lqcd contains two degenerate quark flavours and has a symmetry with respect to the general phase 
rotations in the "two-flavour space": 



<f — > = exp 



(47) 



where to a are arbitrary (x-independent) parameters and a 1 , a 2 , a 3 are the 2x2 Pauli matrices. The 
symmetry transformations (|4*71 form a group SU(2) 5 . 

5 The number of independent parameters for SU(2) is determined in the same way as we did for SU(3) in Lecture 1: one 
counts the number of independent elements in the unitary 2x2 matrix with the unit determinant. 



One has to emphasize that the approximate SU(2)-flavour symmetry emerges "by chance", simply 
because the u and d quark masses turn out to be almost degenerate. Note also that e.m. interaction 
violates this symmetry, due to different electric charges of u and d quarks. But this 0(a em ) effect is 
again small, at the same level of ~ 1%, as the 0( m ^~^ ) violation due to the quark mass difference. 

The approximate degeneracy of u and d flavours manifests itself in hadrons. Replacing u quarks 
by d quarks or vice versa in a given hadron, yields a different hadron which has a very close mass 
and other properties. This qualitative prediction is nicely confirmed by the measured mass differences 
between proton (uud) and neutron (udd), iT + (ud) and ir°([uu— dd]/y2), K + {us) and K°(ds), etc.. The 
typical mass splittings for the u <-> d hadronic partners are at the level of few MeV. Thus, QCD nicely 
explains the origin of isospin symmetry introduced by Heisenberg in 30's in nuclear physics, to describe 
the similarities between the "mirror" isotopes, obtained from each other by interchanging protons and 
neutrons. The second part in the word 'isospin' reflects the analogy with the electron spin symmetry, the 
degeneracy of the spin-up and spin-down electron states in quantum mechanics. 

In fact, one introduces a similar formalism in QCD, attributing isospin 1 = 1/2 to the doublet 
of u and d quarks and treating these two flavours as "up" and "down" components with 1$ = +1/2 
and 7~3 = —1/2, respectively. The hadrons containing u and d quarks in different combinations form 
isomultiplets, with the isospin counting similar to the spin algebra in quantum mechanics. In the case of 
proton and neutron, the diquark ud has isospin 0, therefore, adding u or d quark to the diquark, we get 
the nucleon isodoublet (/ = 1/2) consisting of proton with J3 = +1/2 and neutron with I 3 = —1/2. 
Another isodoublet is formed by K + and K°, where s quark, which has no isospin, is combined with u 
and d, respectively . In the same way, D°(uc) and D~(dc), or B + (ub) and B°(db) build isodoublets. 

Note that antiquarks have the opposite signs of ^3: u (d) has ^3 = — 1/2 (I 3 = +1/2). Combining 
u and d quark with their antiquarks, one gets four states. Three of them belong to isotriplet (/ = 1): 

_ 7/7/ — rjrj 

«d(J 3 = +l), —j=-{h = 0), du{h = -l). (48) 

For example, pions (ir + , it and 7r~), as well as p mesons (p + , p° and p~) form isotriplets. The fourth 
quark-antiquark state is an isosinglet (1 = 0): 



uu + dd 

~7T 



(49) 



which deserves a separate discussion. In general, the uu and dd states transform into each other via 
intermediate gluons. In mesons this transition takes place at long distances, due to some nonperturbative 
mechanism, not necessarily described by diagrams with a fixed number of gluons. In any case, the 
transition amplitude has a characteristic scale of O(Aqcd), much larger than the mass difference ~ 
2(m u — ma) between the uu and dd states. The uu — dd degeneracy yields two orthogonal physical 
states, the ^3 = component of the isotriplet d48t and the isosinglet J49t . Turning to strange quarks one 
encounters the second 1 = state ss. The same gluonic transition mechanism provides a mixing between 
( l49l and the ss-state, Now, the difference between the masses is not small, being of 0(2m s — 2m) ~ 
Kqcd- Hence, the amount of mixing depends on the magnitude of the ss <-> [uu + dd) transition 
amplitude. The latter is quite sensitive to the spin-parity of the meson state. For example, the two 
isosinglet mesons with J p = CP shown in Table |3] have (up to small deviations) the following quark 
content: 

77(547) ~ -^={uu + dd- 2ss) , (50) 
V6 

t/' (958) ~ -^(uu + dd + ss) , (51) 
V3 




Fig. 10: Quark diagram for K ~ > Kir hadronic decays; the four possible combinations of q = u,d and q' = u, d correspond 
to the four decay modes related by isospin symmetry. 



indicating that mixing in the O channel is large. For the J p = 1 mesons the situation is completely 
different. From the quark content of the isosinglet mesons: 

^(782)^^^, 0(1020) -is, (52) 

one concludes that the transition from strange to nonstrange quark pairs in the 1~ state is suppressed. 

Returning to the isospin symmetry, it is worth to mention that it yields useful relations between 
hadronic amplitudes. To give just one simple example, let us consider the four observable hadronic 
decays of K* -meson. The amplitudes of these decays are related via isospin symmetry, so that only one 
amplitude is independent: 

A(K*° -» K+TT-) = — l -A{K*° -> K°ir°) = A(K* + -» K°tt + ) = -^=A(K* + -» K + ir°) . (53) 
V 2 v 2 

To obtain these relations, one does not necessarily need to apply formulae for SU(2) group. In the isospin 
limit, all four decays are described by a single quark diagram, shown in Fig. where the initial and 
final mesons are taken in the valence-quark state. Each individual decay mode has its own combination 
of u and d quarks to be substituted in this diagram. Furthermore, the coefficients 1 / v2 originate from 
the quark content of 7r°. Naturally, the hadronic amplitude attributed to the quark diagram in Fig.lTOlis 
a nonperturbative, long-distance object and cannot be directly calculated in QCD. However all we need 
are the amplitude relations between individual decay modes and not the value of the amplitude itself. 
Squaring the relations (I53t and taking into account the phase space, one predicts the ratios of branching 
fractions confirmed by the experimental values given in [4|. 

2.4 SU(3)-flavour symmetry 

Another flavour symmetry, widely used in hadron phenomenology, is SU(3) n corresponding to the limit 
of QCD with all three quarks, u,d and s, having equal masses. Since in reality m s is only inessentially 
smaller than Aqcd, the magnitude of SU (3) jz-violation in hadrons is not universal, depending on their 
quark content and quantum numbers. 

Neglecting the mass differences between u,d,s quarks and introducing a common mass fh^ one 
approximates Lqqd as 

Lqcd =i L { £= d = s) = * 3 (D^ - m 3 )^3 + L glue + .. , (54) 

where ^3 is a triplet: 




The modified QCD with L^cd ^ has a symmetry with respect to the transitions between 3 flavour 
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(55) 



Although physically, SU(3)ji and the fundamental 5f7(3)-colour have completely different origin, the 
group-theoretical formalism of both symmetries is the same. In particular, the eight A a -matrices entering 
d55t are already given in ©. 

The SU(3)fi symmetry is very helpful in "organizing" the spectra of strange and nonstrange 
hadrons in multiplets. The light-quark meson multiplets are obtained by combining quark flavour-triplets 
and antiquark flavour-antitriplets. Without using the specific rules of SC/ (3)-algebra, which can be found 
in many textbooks, it is easy to figure out that the nine quark-antiquark states split into a singlet ^3 ^3 
and octet ^A ^. The octet has its own isospin substructure 6 . The singlet-octet pattern provides a 
reasonable description for pseudoscalar mesons, in particular, rf meson is close to the SU(3)fi singlet 
state d5TT >. The isotriplet of pions (7r + , it , n~), two isodoublets of kaons (K + , K and K°, K~) and the 
isosinglet 77, given by d50l . together form an octet. However, this pattern is not universal. E.g., in the 
case of vector mesons, uj and <fi states in d52l are neither octets, nor singlets. To complete the counting of 
meson SU(3)fi multiplets, one has to mention also triplets and antitriplets of heavy-light mesons. For 
example, in the case of c quark, D°, D~ , D s (D°, D~ , D s ) form a triplet (antitriplet). 

Generally, SU{3) fi works quite well for baryons, because their characteristic mass scale is a few 
times larger than Aqcd- Importantly, the spin-parity of the lowest baryon SU(3)fi multiplets is fixed, 
due to the total antisymmetry of the baryon "wave-function" in the SU(3) n /spin/colour coordinates 
required by Fermi-statistics. There is an octet with J = 1/2 (including proton and neutron) and decouplet 
with J = 3/2. Let us, for instance, have a look at the latter. It contains the isoquadruplet (/ = 3/2) of 
A-resonances (A ++ (mm), A + (uud), A°(udd), A~(ddd)), the isotriplet of S resonances (E + (uus), 
H°(uds), Y*~(dds)), the isodoublet of H-resonances (E°(uss), E~(dss)) and the isosinglet Q(sss). 
Consulting [4] for the masses of these baryons, one notices a distinct hierarchy: each constituent s quark 
adds an amount of 0(m s ) to the baryon mass. 

Returning to the quark diagram in Fig. \W\ we may now replace the s quark by u or d quark. In 
SU(3) ji limit the new diagram obtained after this replacement is equal to the one with s-quark, yielding 
relations between the K* — ► Kir and p — > mr hadronic amplitudes, e.g. 

A{K* + -» K°7T + ) ~ ^A(p+ -» TT°7T + ) . (56) 

v2 

Such relations are typically violated at the level of 20-30 %, but are still useful from the phenomenolog- 
ical point of view. 



2.5 Heavy quark symmetry 

To complete our survey of flavour symmeties, we now turn to the {c, b} quark sector of QCD Lagrangian. 
The fact that m c ^ S> Aqcd, together with the flavour-independence of quark-gluon interactions, allows 
one to consider an interesting limit of Lqcd where both c and b quarks have infinitely heavy mass: 

m c ~ m b ~ mQ — >■ 00 . (57) 

At first sight, the limit is not justified, because in reality m& is substantially larger than m c . As we shall 
see, the fact that both masses are large turns out to be more important for QCD dynamics. Formally, in 
the limit (B71 one can introduce a doublet of heavy-flavour fields 

»^( S ) ■ • - <*»*> • 

6 In mathematical terms SU(2)isospin is a subgroup of SU(3) 



and rewrite the Lagrangian in a form invariant with respect to SU(2)-rotations in the c, b flavour space: 

Lqcd= J2^Q( iD ^ -rn Q )^ Q + L glue + L uAs . (58) 
Q=c,b 

This particular form of the heavy-quark limit for Lqcd is however not convenient, because the heavy 
mass scale tuq is still present explicitly. To understand why it is desirable to effectively remove that scale, 
let us consider the heavy-quark limit d57l for D or B meson. Since m c and rra& have no direct relation to 
QCD, it makes sense to discuss a hypothetical heavy-light meson H with a mass vtlh containing a heavy 
quark with an arbitrarily large mass uiq. In the rest frame of H the constituent heavy quark stays almost 
at rest, providing a static source of color charge which emits and absorbs gluon fields. The meson mass, 
to a good approximation is 

m H = rrtQ + A , (59) 

where A is the energy of the light quark-gluon "cloud" surrounding the heavy quark. The situation very 
much resembles the hydrogen atom where the total mass of the atom is a sum of an extremely large (~ 
GeV) proton mass and a small energy of the electron cloud (~ MeV). The essential point is that the 
electron itself is nonrelativistic. One can isolate the electron mass from the rest of the energy, introduce 
the Coulomb potential and kinetic energy, and eventually solve the equations of motion, determining the 
electron energy levels. In heavy hadrons the light-quark cloud is purely relativistic {m u ^, s < ^-qcd) 
and has a complicated long-distance nature. 

Nevertheless, one essential feature is common for both bound states. In the atom the energy of the 
electron cloud does not depend on the proton mass. Likewise, in the heavy-light meson A in d59l is (up 
to 1/niQ corrections) independent of rriQ. From atomic physics we know that the electron energy levels 
in hydrogen and deuterium coincide to a great precision. The fact that deuteron is twice more massive 
than proton, does not play a role for the energy levels, because in both cases the atomic nuclei are static. 
Important is that the electric charge does not change by switching from proton to deuteron. Similarly, in 
the H meson A changes very little if one replaces uiq by or by m c , because the colour charge of the 
heavy quark does not change. Thus, the heavy-flavour symmetry is in reality the symmetry between the 
light-quark remnants of the heavy hadrons, so that the heavy-quark mass scale indeed plays a secondary 
role. 

To achieve a quantitative level, a special formalism of heavy quark effective theory (HQET) was 
developed for applications of QCD to heavy-light hadrons. One starts from Lqcd and introduces trans- 
formations which decouple the ~ mg part of the heavy-quark field from the part which has the remnant 
momentum ~ Kqcd- Only the latter part strongly interacting with the light quark-gluon cloud is relevant 
for QCD dynamics. In HQET one integrates out the heavy degrees of freedom and works with the La- 
grangian containing a new effective quark field carrying the flavour of Q but no mass. Not only the static 
limit d57b but also an expansion in powers of 1/ iuq can be systematically treated. The field-theoretical 
aspects of heavy-mass expansion are nicely explained in the literature (see e.g., [8], [9]); I will only focus 
on some important phenomenological consequences of heavy-flavour symmetry. 

One famous example is the B — > Dlv\ decay involving weak b — > c transition in Standard Model. 
The unknown part of the decay amplitude is the hadronic matrix element 

(D(p D )\c^b\B(p B ) (60) 

determined by the long-distance interactions involving the initial and final heavy quarks as well as the 
surrounding light quark-gluon "cloud". One chooses a special kinematical configuration, the "zero recoil 
point" where the momentum transfer to the lepton pair is equal to 



q 2 = (PB - Pd) 2 = {m B - m D ) 2 



(61) 



In the B meson rest system p B = (m#, 0, 0, 0) this point correspond to the final D meson at rest. In the 
heavy-quark limit the replacement of b quark by c-quark does not change the hadronic state: 

(D(p D )\ = (B(p B )\ 

and the matrix element (161 b reduces to a trivial normalization factor. One can therefore predict the decay 
amplitude in the zero recoil point up to 1/mg corrections. In fact there is a theorem stating that these 
corrections are even smaller and start from 0(l/mg), but to derive this and other important details one 
needs a full-scale HQET framework. 

Without resorting to the effective theory, it is possible to understand the origin of another important 
symmetry emerging in the heavy-quark limit. In the hydrogen atom, the electron and proton have mag- 
netic moments related to their spins and yielding interactions with the external magnetic fields or with 
each other (spin-spin interactions). The magnetic moments are inversely proportional to the masses, so 
that the proton magnetic moment plays no role for the electron energy levels. Each level is degenerate 
with respect to the proton spin direction. Since QED and QCD have very similar vector boson inter- 
actions with spin 1/2 particles, the spin 1/2 quarks also have chromomagnetic moments, and interact 
with the "magnetic" parts of gluonic fields and with other quarks. For the heavy nonrelativistic quark 
the chromomagnetic moment is inversely proportional to the heavy quark mass ttiq. In the infinite mass 
limit the interaction vanishes and, hence the light-cloud energy A is independent of the spin orientation 
of the heavy quark. 

One arrives at a new classification of heavy-light states based on this heavy-quark spin symmetry. 
Instead of adding together the orbital momentum and the total spin of quarks as we did in d42l it is more 
appropriate to introduce, for a Qq meson (Q = c,b; q = u, d, s), the total angular momentum of light 
degrees of freedom: 

Jlight = L + S q (62) 

Adding the heavy quark spin sq = 1/2 to Ju g ht one gets degenerate doublets of heavy-light mesons 
with total angular momentum J = Ju g ht ±1/2. At L = one simply has Ju g ht = 1/2 and therefore 
a doublet of mesons with J p = 1~ and J p = 0~ consisting of B and B* (D and D*) in the b quark 
(charm) sector. The mass differences within doublets are indeed very small |4) : 

S B = m B * - m B = 47 MeV, 5 D = m D * - m D = 142 MeV, (63) 

indicating that the heavy-quark spin symmetry works quite well, especially for the heavier b quark. 
Taking into account that the mass differences are ~ l/ m Q effects, one expects that 5 b /5d — {m c /mb) 
which is also in accordance with ( I63t and d3"Tl) . I leave as an exercise to show that at L = 1 there are two 
degenerate doublets: one with J p = + , 1 + and another one with J p = 2 + , 1 + . 

2.6 Exotic hadrons 

The colour confinement principle does not exclude hadronic states with an "exotic" valence quark con- 
tent, different from qq or qqq. Quarks, antiquarks and gluons can be added together in any combination, 
e.g., qqG , GG , qqqq or qqqqq, provided they are in a colour-neutral state. Since one cannot calculate 
the spectrum of hadrons in QCD with a good precision, predictions of exotic states are generally model- 
dependent. It is always problematic to distinguish an exotic hadron from the excitation of an ordinary 
hadron with the same J p and flavour quantum numbers. Moreover, in this case one expects mixing 
between ordinary and exotic hadrons. For example, if there is a J p = 0~ state composed of two gluons 
GG (glueball), it should be mixed with rj to a certain degree, so that r( acquires a glueball component. 

Therefore, the most interesting, "smoking gun" signatures are the hadrons with exotic quantum 
numbers (flavour content and/or J PC ), forbidden for quark-antiquark mesons or three-quark baryons. 
For example it is impossible to arrange a flavour-neutral quark-antiquark state with J PC = 1 h . The 
P and C ( charge-conjugation) parities of a fermion-antifermion state are determined by the rules: P = 



Fig. 11 



A schematic view of pion-proton elastic scattering at large momentum transfers. 



— (— 1) L and C = (— 1) L+S , so that P = — 1 means L = 0, 2, .... Hence, the only possibility to have 
C = +1 is S = 0, 2, ... However, adding together even values of L and S one cannot get J = 1. On 
the other hand, adding one constituent gluon to a qq pair, one easily makes a "hybrid" meson with qqG 
content and J PC = 1 quantum numbers. Searches for hybrid mesons are currently being carried out, 
but the experimental situation is not yet settled. 

The recently observed narrow baryon resonance 0(1540) decaying to K + n\ 10] is another promis- 
ing candidate for hadron exotics, a state suudd with five valence constituents (pentaquark). Flavour 
symmetries are important model-independent tools to confirm/reject the experimental candidates for ex- 
otic resonances. In particular, an important task is to find the symmetry partners of these hadrons and 
to fill the relevant isospin and SU(3)ji multiplets (in case of the pentaquark it is actually the SU(3)fi 
antidecouplet). 

3. QCD AT SHORT DISTANCES 

3.1 Probing short distances with electroweak quark currents 

In this lecture we return to the quark-gluon interactions at large momentum transfers (short distances). 
In this region, practically at Q > 1 GeV, a s is small, allowing one to apply the perturbative expansion. It 
is then possible to test QCD quantitatively, calculating various quark-gluon interaction processes at large 
Q and comparing the results with the available experimental data. Note however, that the tiaditional way 
to study interactions by scattering one object on the other is not applicable to quarks and gluons. They 
simply are not available in free-particle states. One needs to trace quarks inside hadrons, where the long- 
distance forces are important. Take as an example the elastic pion-proton scattering at large momentum 
transfers (see Fig. fTTTi . Here one has to combine the perturbative quark-quark scattering amplitudes 
(two-gluon exchange) with the "wave functions" of quarks inside the initial and final hadrons. To obtain 
these functions one needs to go beyond perturbative QCD. Therefore, an unambiguous extraction of the 
perturbative amplitude from the data on the scattering cross section is not a realistic task. 

The situation is not so hopeless, actually, since we have at our disposal electroweak bosons 
(7, W, Z) interacting with quarks in a pointlike way. Electroweak interactions at large momentum trans- 
fer serve as external probes of short-distance dynamics. In Lecture 1, we already discussed one example: 
the quark-antiquark pair production in Z decay. To list all possible electroweak sources of quarks in a 
more systematic way, I start with the photon. Its interaction with the quark e.m. current was already 
given in ©, let me write it down again: 




(64) 




introducing a compact notation for the quark e.m. current: 



(65) 



q=u,d,s,c,. 



where the summation over colour indices is not shown for brevity. The quark weak current entering the 
quark- W flavour-changing interaction: 



L w (x) 



is more complicated and includes the CKM mixing matrix: 
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Finally, the quark-Z interaction is 



L z 



9 -z 7 
2cos9 w ^ ' 



where the quark electroweak neutral current is a mixture of vector and axial-vector parts. 



(66) 



(67) 



(68) 



3.2 Perturbative QCD in e + e — > hadrons 

In e + e~ annihilation at high energies the virtual photon provides a short-distance source of quark- 
antiquark pairs. This process is hadron-free in the initial state. The photon-exchange diagram is depicted 
in Fig. [^] (for simplicity, I ignore the additional Z-exchange diagram). The experimentally measured 
total cross section <j to f(e + e - — > h) depends on one kinematical variable s = (p e _ + p e +) 2 - The virtual 
timelike photon transfers its energy to the hadronic state. At very large y/s 3> Aqcd and yfs S> m q , 
the initial pair of quarks is produced at an average distance of 0(l/^fs), much smaller than the typical 
hadronic distance scale 1/Kqcd- Due to asymptotic freedom of QCD, gluonic interactions of the pro- 
duced quark pair are suppressed by small a s . Hence, e + e~ annihilation at high energies provides an 
almost pointlike source of quasi-free quark pairs. At long distances the created quarks and antiquarks 
are inevitably converted into some hadronic state. Since in the total cross section the summation is done 
over all hadronic states produced at a given energy, the total probability of hadronization sums up to a 
unit. Hence, at y/s — ► 00 the hadronic cross section is well approximated by the cross section of the free 
quark-antiquark pair production: 

4i e '" h \s)^ E * (0+9 "-^W, (69) 

q=u,d,s,.. 



summed over all quark flavours with m q <C \fs. This, so called parton model approximation for e + e~ — > 
h is confirmed by experimental data. Moreover, the majority of events saturating the cross section at high 
yfs consists of two distinct hadronic jets originating from the initial quark pair. 

The way we obtained (15^1 may seem too qualitative and a bit "hand-waving". In the following, we 
shall derive the asymptotic cross section ( 169 \ in a more rigorous way. In this derivation several important 
concepts will be introduced, to be used in discussing further topics covered by these lectures. 

We start with the formal definition of the total cross section: 



( e + e _>/!) _ 
hot W ~ 2g 



' ' J2\( h n\T\e + e-)\ 2 , (70) 



where the sum over the final hadronic states h n includes phase-space integration and implies summation 
over spins (polarizations). The matrix elements (f\T\i) = Tfi (in (f70t \i) = \e + e~) and (/| = {h n \) 
determine the general S'-matrix of the theory: 

S fi = (f\S\i)=S fi + iT fi . (71) 

The usual representation of S'-matrix in terms of Lagrangian has the time-ordered exponential form: 



S = rjexp i J d i x{L QCD {x) + L QED (x)) 



(72) 



where Lqed includes e.m. interactions of quarks and leptons. Furthermore, the unitarity of S matrix is 
used: 

SS f = 1 , (73) 



or 



J2(f\S\n)(n\S^\i)=6 fi . (74) 



From now on we consider the forward scattering f = i. Replacing (n\S^\i) = (i\S\n)* and substituting 
dm in d74l . one obtains the unitarity relation for Ta (the optical theorem): 

2\mT u = ^ |T ni | 2 . (75) 

n 

To apply this universal relation to e + e~- scattering, we take |i) = \e + e~) with four-momentum q = 
Pe+ + Pe-> so tnat q 2 = s > 0. Furthermore, we choose |n) = \h n ) restricting the set of intermediate 
states by hadronic states. As a result we obtain a rigorous relation between the amplitude of the forward 
e + e~ — * h — > e + e~ scattering via hadronic intermediate states, Ta = A^ e+e -* h ~* e+e ) ; and the sum 
over the squared e + e~ — > h n amplitudes. The latter sum, according to (f70t is proportional to the total 
hadronic cross section. The optical theorem d75t takes the form: 

2Im^ e+e -^ e+e -)( S ) = ^|(/ ln |T|e + e-)| 2 = 2 S a^"^ ) ( S ). (76) 

Diagrammatically, this relation is represented in Fig. [O] The amplitude 

A(e+e -^ e+e - ){g2) = _^_ ( ^ 7 ^ e)( ^ 7>e)n ^ (g) (77) 

contains the photon propagators and the products of electron and positron spinors in both initial and final 
states, written according to the QED Feynman rules. The nontrivial part of this amplitude, denoted as 
li^y, describes the — ► h — ► jl m transition, and is called the correlation function (or correlator) of 
quark currents. The formal expression for this object reads: 

IVfa) = i f d 4 x e^(0\T{jf l m (x)jt m (x)}\0) . (78) 




Fig. 13: Unitarity relation for e + e —> hadrons 



Due to the conservation of e.m. current (d^j^ 71 = 0), the correlation function depends on one invariant 
amplitude: 

iV(<z) = (-M 2 + ^)n(<? 2 ) . (79) 

Substituting (1791 in (1771 and taking imaginary part from both sides, we obtain 

Im^ e+e ^^ e+e_ )( S ) = --(^^e)(^e7^ e )Imn( S ) . (80) 

s 

It is convenient to normalize the hadronic cross section to the e + e~ — > fL + fi" cross section known 
from QED. One can literally repeat the derivation done above, taking instead of hadronic states the 
states: \n) = The resulting relations are quite similar to dToT ). dTTl and dSDl : 

2Im^ e+e "^ + ^"^ e+e ")(s) = Y, \(v + V~\T\e + e-}\ 2 = 2sa (e+e "^ + ^" ) (s) , (81) 

<A (e+e-- + M + A.--« + e-)( g 2) = i-(^ e )(^ e )Il0j)( ? ) , (82) 

and 

Im^ e+e "^ + ^"^ e+e -)( S ) = ~$ e We)$elM*^\8) , (83) 

where the muonic correlation function is nothing but a 2-point muon-loop diagram. Furthermore, 
the cross section in (f8TT> taken from QED textbooks, reads: 

(j(e + e -_ >M + M -)( s ) = i^lm . (84) 
3s 

Dividing the hadronic unitarity relation dTBT i by the muonic one (ISTl and using dSOl and ( l8*3l . we obtain 
a useful ratio: 

Imn( S ) 



ImnW(s) ^(e+e-^+M )(f 



i?(s) . (85) 



The next key point in our derivation is the analysis of the correlation function II (q 2 ) at spacelike 
q 2 < 0. At large \q 2 \ = Q 2 » A.qcd< tne long-distance domain in the space-time integral in ( TTSl is 
suppressed by the strongly oscillating exponent and the short distances/times |x| ~ xo ~ 1/Q dom- 
inate. This justifies using QCD perturbation theory with a s (Q) —> at Q — > oo. The leading-order 



(a) 




(b) (c) (d) 

Fig. 14: Diagrams corresponding to the perturbative contributions to the correlation function of two quark currents: (a) the 
leading-order loop and (b-d) the 0(a s ) corrections. Wavy lines denote external currents with 4-momentum q, solid lines quarks 
and dashed lines gluons. 



asymptotically free result is given by the 2-point quark-loop diagram (Fig.[T3h). and the next-to-leading 
corrections are determined by 0(a s ) two-loop diagrams (Fig. [T4b.c.d). Calculation of these diagrams 
yields (at m q <C Q): 




Ii(q 2 ) = -—r \ Qi log Mr U + — )+° " )+ const, (86) 



n (5 2 ) = — A I . (87) 



where each flavour q contributes with the same expression and a coefficient Q 2 . The natural scale for a s 
is q 2 . In recent years, due to tremendous calculational efforts, the 0{a 2 ) and even O(a^) corrections to 
n have been calculated; also the loop diagrams for the massive quark are known with a high accuracy. 

In the final stage of our derivation the calculated function n(g 2 ) at q 2 < is related to Imn(s) 
at positive s. One employs Cauchy theorem for the function H(z) obtained from H(q 2 ) by analytically 
continuing the real variable q 2 to the complex values, q 2 — > z: 

h 

2m J z — q z 

c 

The integration contour C is shown in Fig.^] It circumvents the singularities of the function 11(2;), i.e. 
the points/regions where Imn(z) / 0. According to dTBT ) and dSUl the location of singularities at real 
q 2 > is determined by the masses of resonances and/or the thresholds of multiparticle hadronic states 
produced in e + e~ — > h; the lowest one is at s m i n = Am 2 corresponding to the threshold of the lightest 
two-pion state. Subdividing the contour C into: 1) a large circle with the radius R, 2) an infinitely small 
semicircle C surrounding s m i n and 3) two straight lines from s m i n to R, we can rewrite the integral in 
terms of three separate contributions: 

„, 2s 1 f , n(z) 1 f n(z + id) -IL(z-iS) 1 t n(z) 

|z|=_R *min Q 

Suppose the function decreases at \q 2 \ — ► oo: H(q 2 ) ~ l/|g 2 | A , where A > 0. Then the first integral 
vanishes at R — > oo. Taking infinitely small semicircle, one makes the third integral also vanishing. 
Furthermore, since there are no singularities of n(z) at Rez < s m i n , the integrand in the second integral 
reduces to the imaginary part: H(q 2 + i5) — H(q 2 — i5) = 2ilmll(g 2 ) (due to Schwartz reflection 



Fig. 15: The integration contour in \87\ . The crosses on the positive real axis indicate singularities ofH(q 2 ) 



principle). Finally, we obtain the desired dispersion relation 



oo 

o If, ImU(s) 



Il(q z ) = - / ds (89) 

7T J s — q z — id 

with l.h.s. calculated in QCD in a form of perturbative expansion d8*6l l and r.h.s. related to the cross 
section via (f85T >: 

ImII(s) = R(s)lmU M (s) . (90) 

The quantity R(s) is directly measurable in e + e~ experiments. It remains to determine Imll^^s). 
Knowing the answer for the quark loop diagram in Fig. IT4lt it is very easy to write down the expression 
for the muon loop. Since we neglect masses in both diagrams, the only difference is the factor 3, from 
summing up the colour states in the quark loop. This factor is naturally absent for the muon loop. From 
(186b . taking imaginary part, one obtains 

ImnW( S ) = -^-. (91) 

1Z7T 

Finally, to guarantee the convergence of the dispersion integral d89l ). 7 let us differentiate both parts in 

I 2 : 



dq 2 ir J (s - "'- v ~ 

where d86t gives for l.h.s. 



9 



Note that divergent and constant terms disappeared after differentiation and play no role in our derivation. 
The final form of the dispersion relation is 



it J J [s — q 



1 + -+...) = / d ST ^L-, (94) 



7 For brevity, I avoid a longer derivation which includes some special mathematical construction (subtractions). 




where ellipses denote higher-order in a s corrections. The fact that d94b is valid at (— q 2 ) — > oo unam- 
biguously fixes the constant limit of R(s) at large s: 

n f 

R{s)^3^Q 2 q , (95) 
q 

where nf indicates that R includes all "active" quark flavours, for which the condition S> m q is 
fulfilled. Finally, we notice that (l95l coincides with the parton model prediction d69b . taking into account 
that the free-quark- and muon-pair cross sections differ only by the colour factor times the quark charge 
squared: 

Importantly, QCD not only reproduces the parton model prediction for R(s) but also provides perturba- 
tive corrections, as well as predicts the integral d94b over R(s). 

The experimental data collected in various regions of \/~s nicely confirm d95t . According to FiglTol 
taken from [4], the ratio R(s) approaches first the constant value R ~ 2 at energies ~ 2 — 3 GeV, above 
the region of vector meson resonances p,ui,(j) (and below the charm threshold). That is exactly the value 
anticipated from d95b for nj = 3. Well above charmonium resonances, a new constant level is achieved: 
R = 2 + 3Q 2 = 10/3. And finally, R = 10/3 + 3Q 2 , = 11/3 is settled at energies above T resonances, 
where all five quark flavours are in their asymptotic regime. Actually, the current data on R(s) are so 
precise that one should also include small a s -corrections to Rqcd- 

There are other similar inclusive observables calculable in QCD, among them the total widths 
^tot{Z — ► hadrons) and T to t (W — ► hadrons). They have the same status as R(s), but a fixed scale mz 
or mw instead of y/s. One has also to mention an interesting and well developed sub-field of perturbative 
QCD related to the jet and/or heavy-quark production in e + e _ and hadron collisions at high energies. The 
underlying short-distance quark-gluon processes are successfully traced in the experimentally observed 
multijet structure of the final state. Naturally, hadrons cannot be completely avoided, because, after 
all, quarks and gluons hadronize. In fact, hadronization in jet physics is nowadays a somewhat routine 
procedure described by QCD-oriented models (e.g. the Lund model integrated within PYTHIA [11]). At 
lower scales, Q ~ 1 — 2 GeV, inclusive decays of r-lepton are among useful tools to study perturbative 
QCD (see, e.g. El). 



3.3 Deep inelastic scattering and operator-product expansion 

The processes of lepton-nucleon deep inelastic scattering (DIS), IN — > Ih and viN — > Ih (I = e, a), 
are among the best testing grounds of QCD. The long-distance structure of the initial nucleon makes 
these processes more complicated than e + e~ — > hadrons. For definiteness, let us consider the electron- 
nucleon scattering mediated by the virtual photon 7* with the 4-momentum q, whereas p is the nucleon 
4-momentum. To measure the DIS cross section, one only has to detect the final electron. In the nucleon 
rest frame the invariant variables Q 2 = —q 2 and v = q ■ p are related to the initial and final energies and 
the scattering angle of the electron: Q 2 = AEE'sin 2 6, v = (E — E')m^. 

The specific kinematical region Q 2 ~ v S> A.q CD has to be chosen to reveal the spectacular effect 
of asymptotic freedom. In this region the experimentally measured DIS cross section, normalized to the 
cross section of the electron scattering on a free pointlike quark a po i n t, depends essentially on the ratio 
x = Q 2 /2v. This effect (Bjorken scaling), was first interpreted in terms of a beautiful phenomenological 
model suggested by Feynman. One considers the reference frame with large nucleon momentum p, so 
that — q gives the longitudinal direction. Neglecting long-distance binding forces between quarks, the 
initial nucleon is represented as a bunch of free constituents {partons): quarks, antiquarks and gluons 
moving in the longitudinal direction 8 . For simplicity the quark masses and transverse momenta as well 
as the nucleon mass are neglected, in comparison with Q 2 and v. The electron scatters on one of the 
nucleon constituents (excluding gluons, of course, because they are electrically neutral) which has the 
momentum fraction £ , so that after the pointlike collision the quark 4-momentum is p^ + q. The massless 
quark has to remain on-shell, (£p + q) 2 = 0, therefore 2£(p ■ q) + q 2 = and 

£ = Q 2 /2v = x . (96) 

The cross section is then represented as a sum of all possible "elementary" processes integrated over £ 

1 

da (eN^eh) {Q 2^ v) = J ^ f^) da^Q 2 , v)5{£ - *R = £ da^Q 2 , u)f t (x) , (97) 

i=u,u,d,d... i 

where 5-function takes into account d96l l. The sum goes over all quark and antiquark species inside 
proton, and /«(£) are the partem distributions defined as the probabilities to find the z-th constituent with 
the momentum fraction £. 

In QCD, the scaling behavoiur d97t corresponds to the asymptotic-freedom approximation (a s = 
0). Switching on the perturbative quark-gluon interactions one predicts logarithmic corrections ~ log(Q 2 ) 
to this formula which are nicely reproduced by experiment. There are many excellent reviews and lec- 
tures where DIS in perturbative QCD are discussed (see e.g.[13|). Here we shall concentrate on one 
essential aspect: separation of long- and short-distance effects. 

Let us approach DIS from the quantum-field theory side, as we did above for e + e~ — > hadrons. 
Omitting for simplicity the initial and final electrons, one can represent the hadronic part of the DIS cross 
section in a form of 7*iV — > h cross section. Employing unitarity relation d75t with \i) = |7*iV) one is 
able to relate the DIS cross section to the r y*(q)N(p) — ► r y*(q)N(p) forward-scattering amplitude (see 
Fig. EH: 

2lmA^ N ^ h ^ N \q,p) =Y,\(hn\T\j*N}\ 2 ~ do^* N ^ h \Q 2 ,v) . (98) 

The amplitude 

A^ N - h ^* N \q,p) = e»(qy(q)T^(q,p) (99) 
contains the photon polarization vectors multiplied by a new purely hadronic object 

2^(p, q)=ij d 4 xe^(N(p)\T{jfr(.x)jZ m (0)}\N(p)} , (100) 

8 We refer here to a generic picture of the nucleon, where all possible multiparticle components are coherently added to the 
valence three-quark state, similar to <4U for the pion. 




Fig. 17: Unitarity relation for DIS cross section 



which resembles the correlation function of two currents we introduced above, but, instead of vacuum, 
has nucleons in the initial and final states. 

One can prove that at large Q 2 and v the dominant contribution to the above integral stems from 
small space-time intervals x 2 ~ l/Q 2 ~ l/v 9 . In other words, the points of absorption and emission of 
the virtual photon are located close to the light-cone x 2 = 0. The process takes place in the asymptotic 
freedom regime, that is, a single quark absorbs the photon, and penetrates quasi-freely at small x 2 before 
emitting the photon. All other contributions, for example, with different quarks emitting and absorbing 
initial and final photons are suppressed by powers of l/Q 2 , l/v. To describe the free-quark propagation 
from x to we use the quark propagator 

(0|T{^(x)^(0)}|0) = J dtp (^f) e-*« = • (101) 

neglecting the quark mass. Substituting the e.m. currents in d99l in terms of quark fields, contracting the 
fields of the propagating quark and using dlOlt we obtain 

T^(p,q)=ijd 4 xe^ Yl ( N (p)\nMxhM x )^hu^ q (0)}\N(p)) 

q=u,d,s,.. 

= i J d4xeigX ^7^2 E <iV(p)|^(x) 7M 7a7^,(0)|JV(p)> + ... (102) 

\ ' q=u,d,s,.. 

In this expression, where only the leading term is shown, the calculable short-distance part (the quark 
propagator) is separated from the long-distance part which is represented by the quark-antiquark matrix 
element taken between nucleon states. This hadronic matrix element is a complicated object which has 
to be either determined from experiment or calculated using methods beyond perturbative QCD which 
will be discussed in the next two lectures. I skip the derivation of the cross section formula from the 
imaginary part of T^ v which allows to relate the matrix element introduced in dl02t with the parton dis- 
tributions. Also a s corrections can be systematically calculated; they contain important and observable 
logQ 2 effects. Important is that the long-distance matrix element (or parton distributions) are universal 
characteristics of the nucleon and they do not change if the short-distance part of the process changes 
(e.g., by switching to neutrino-nucleon scattering W*N — ► h c where the charmed quark is produced). To 

'Note that the light-cone dominance x 2 — x 2 t — x 2 ~ is a more general condition than the small-distance/time dominance 
xq ~ x ~ which takes place in e + e~ —> hadrons. 



summarize, in DIS it is possible to separate the short-distance domain by choosing appropriate kinemat- 
ical region and defining convenient physical observables. The short-distance quark-gluon interactions 
are calculable in a form of perturbation theory in a s , whereas the long-distance part is parameterized in 
terms of universal hadronic matrix elements. The procedure of approximating the product of currents in 
(I102t by a quark-antiquark operator and separating short and long distances is called operator-product 
expansion (OPE). Implicitly or explicitly, OPE is used in almost any perturbative QCD treatment of 
hadronic processes. 



4. LONG-DISTANCE DYNAMICS AND QCD VACUUM 
4.1 Vacuum condensates 

QCD at short distances does not essentially help in understanding the long-distance dynamics of quarks 
and gluons. From the short-distance side we only know that the running coupling a s (Q) increases at 
low-momentum scales and eventually diverges at Q ~ Aqcd (see Fig. 0. Is the growth of a s the 
only dynamical reason for confinement? It is not possible to answer this question remaining within the 
perturbative QCD framework, because the language of Feynman diagrams with propagators and vertices 
is not applicable already at a s ~ 1. QCD in nonperturbative regime is currently being studied using 
other methods, first of all, the lattice simulations. From these studies, there is a growing confidence that 
long-distance dynamics is closely related to the nontrivial properties of the physical vacuum in QCD. 

For a given dynamical system, vacuum is a state with the minimal possible energy. Evidently, in 
QCD the vacuum state contains no hadrons, because creating any hadron always costs a certain amount of 
energy. In that sense, QCD vacuum has to be identified with the (0 1 state in the correlation function (I78t . 
Given that the vacuum state contains no hadrons does not yet mean that it is completely empty. There 
could be quantum fluctuations of quark and gluon fields with nonvanishing densities. The existence of 
vacuum fields is manifested, e.g. by instantons, special solutions of QCD equations of motion having a 
form of localized dense gluonic fields (for an introductionary review on instantons see, e.g., 031)- Lattice 
QCD provides an independent evidence for quark/gluon fields in the vacuum. Without going into further 
theoretical details, I will rather concentrate on the phenomenology of vacuum fields in QCD. We shall 
see that properties of hadrons are influenced by the existence of quark and gluon vacuum fluctuations 
with nonvanishing average densities, the so called vacuum condensates. 

Formally, in the presence of vacuum fields, the matrix elements of quark and gluon field operators 
between the initial 1 0) and final (0 1 states are different from zero. The combinations of fields have to 
obey Lorentz-invariance, colour gauge symmetry and flavour conservation, so that the simplest allowed 
composite operators are 

3 = Vv?/4, 4 = G%G^\ 5 = ^ qi (X a )la^G a ^^ k q , 



Oa 



where = (i/2) [7^, 7^] and T a are various combinations of Lorentz- and colour matrices. The indices 
at Od reflect their dimension d in GeV units. Furthermore, if the operators are taken at different 4-points, 
care should be taken of the local gauge invariance. For instance, the quark-antiquark nonlocal matrix 
element has the following form: 

(0| %(x)[x,0]il> q (0) |0>, (104) 

where [x, 0] = exp ig s Jq dvx' 1 A'^(vx)(X a /2) is the so-called gauge factor. Only the matrix elements 
with the light quarks q = u, d, s are relevant for the nonperturbative long-distance dynamics. A pair of 
heavy c (6) quarks can be created in vacuum only at short distances/times of 0(l/2m c ) (0(l/2mi,)), 
that is, perturbatively. 

Without fully solving QCD, very little could be said about vacuum fields, in particular about their 
fluctuations at long distances which have typical "wavelengths" of 0(1/ 'Aqcd)- Thus, we are not able 



to calculate the matrix element (11041) explicitly, as a function of x. It is still possible to investigate the 
vacuum phenomena in QCD applying certain approximations. One possibility is to study the average 
local densities. The vacuum average of the product of quark and antiquark fields, 

<0|^*|0> = (gg>, (105) 

corresponds to the x — ► limit of the matrix element (U04t . The simplest vacuum density of gluon fields 
is 

(0 | G a ^G a ^ | 0) = (GG) . (106) 

Due to translational invariance, both (qq) and (GG) are independent of the 4-coordinate. These universal 
parameters are usually called the densities of quark and gluon condensates, respectively. As we shall 
see in the following subsection, the nonvanishing quark condensate drastically influences the symmetry 
properties of QCD. 



4.2 Chiral symmetry and its violation in QCD 

Let us return to the isospin symmetry limit d46l > of the QCD Lagrangian: 

L QCD = ~*( iD ^ ~ rh)V + Lglue + - • (107) 

Since rh ~ m u ~ <C A.qcd, a reasonable approximation is to put rh — > 0, so that u- and (i-quark 
components of the ^-doublet become massless. 

Each Dirac spinor can be decomposed into the left- and right-handed components, 

I + 75 1 — 75 
*Pq = 5 ^ q H 2 ^ = ^1 R + ^1 L ' ( 108 ) 

where, by definition, the left-handed (right-handed) quark has an antiparallel (parallel) spin projection 
on its 3-momentum. Similarly, for the conjugated fields one has: 

i>q = i>q l +^g 1 ^ 75 =4>qR + i>qL ■ (109) 

Rewriting ^ in terms of the left- and right-handed components we obtain the following decomposition 
of the Lagrangian dl07t in the massless limit: 

l qcd = *R iDn-fVR + *l iD^ L + L glue + ... . (1 10) 

The quark-gluon interaction term in Lqqq is now split into two parts, g s ^>ji'y fl A IJ ' a (X a /2) 1 ^ji and 
9s^L"f[ J ,A lia (\ a /2)*i, so that quarks conserve their chirality (left- or right-handedness) after emitting 
/absorbing an arbitrary number of gluons. Importantly, also the interaction of quarks with photons has 
the same property (take the e.m. interaction d64l and decompose it into L and R parts in the same 
way as above). In the massless (chiral) limit, quarks of left- and right chiralities propagate and inter- 
act independently from each other. In fact, it is possible to introduce two independent isospin SU(2) 
transformations, separately for L and R fields: 



^ L -> V' L = exp (-*y<^) , 



* R - <t>' R = ex V \-i—ui\y R . (Ill) 

Restoring the mass in leads to the violation of chiral symmetry. The Lagrangian mass term can 

be represented as an effective transition between left- and right-handed quarks: 

= m ( ^l^r + ^r^l) ■ (112) 




Fig. 18: One of the diagrams describing the photon-proton scattering. The second diagram is obtained by interchanging the 
photon lines. 



Having in mind the smallness of the it, cf-quark masses, one naturally expects that an approximate chiral 
symmetry manifests itself in the observable properties of hadrons. In reality, the symmetry is violated 
quite substantially, as the two following examples demonstrate. 

The first one is a real hadronic process: the photon scattering on a polarized (e.g. left-handed) 
proton. One of the diagrams is shown in Fig. (fT8l . This process is a very complicated mixture of quark- 
photon and quark-gluon interactions at long distances, determined by the quark structure of the proton. 
Importantly, all these interactions obey an approximate chiral symmetry at the Lagrangian level. Hence, 
the amplitude of the proton-chirality flip is expected to be very small, 0(rn Ui d/ A-qcd) ~ 1% of the 
total scattering amplitude. To check this conjecture, let me consider the case when the initial photon 
energy is much smaller than the proton mass. In this case, it is possible to approximate the proton with 
a point-like particle. The 7p — > 7p amplitude is then simply obtained from the — > amplitude in 
QED (Compton effect), replacing the electron spinors and propagators by the proton ones: 



A(jp — > jp) ~ e 2 u{p2) 



(pi + ^l)^ + m p (pi - k 2 )^ + m 



^ a (pi + ki) 2 - ^ f3 + ^ 13 ( pi _ fc 2 y _ m 2 



UL{pi)<%£, (H3) 



where are the polarization vectors of the initial and final photons with the 4-momenta k\ and k 2 , 
respectively. Without even completing the calculation, we notice that the ~ m p terms in the amplitude 
flip the initial proton chirality L — > R, whereas the ~ 7^ terms preserve chirality. Importantly, the 
contributions of both types are of the same order, determined by the scale m p , indicating that chiral 
symmetry for the photon-proton scattering is broken at 100% level. 

To present the second example of the chiral symmetry violation, I start from the correlation func- 
tion 

U^=ij d 4 xe^(0|T{^(x)>(0)}|0) = (- M 2 + q^uMq 2 ) , (114) 
very similar to the one introduced in Lecture 3, but containing a slightly different quark current: 

U = -A^ulti^u ~ i>dl^d) , (115) 

which produces the 1=1 and J p = 1~ quark-antiquark states. Note that this current is conserved, 
dftjV = 0, even if m Uj d / 0. We then follow the same derivation as in Lecture 3 and obtain the 



dispersion relation (189b for H(q 2 ). The only change is in the imaginary part d90l . where now only the 
states with 1 = 1 contribute to the total cross section (p meson and its radial excitations, the two-pion 
state and other states with an even number of pions), so that the ratio R(s) has to be replaced by 



a{e + e — ► fi + p ) 

In parallel, we consider the correlation function 

11%, = ij d 4 ^(0|T{i M5 (x)> 5 (0)}|0) (117) 
of two axial- vector currents with 1=1: 

This current is conserved only in the chiral symmetry (m u ^ = 0) limit: 

= — ^(2m«V)„75'0 u - 2m d i>dlhi'd) ■ (119) 
Decomposing the correlation function (11171) in two tensor structures: 

njU?) = ~9^{q 2 ) + wJktf) > ( 12 °) 

we notice that in the chiral limit there is only one independent invariant amplitude: 

n 5 (g 2 ) = ^. (121) 



The dispersion relation for II5 has the same form as ( 1891 ): 

00 

n 5(? 3 ) = I Jds^L, (in) 

7T J s — q z — 10 



s 5 . 

mi a 



where is the corresponding threshold. In order to determine the imaginary part via relation similar to 
d90h we introduce a slightly artificial cross section <7 to4 (e + e~ — > Z A — > h(I = 1)) of hadron production 
mediated by Z boson coupled to the axial-vector current. Then, 



lmU 5 (s) = -^R { 5 I - l) (s), (123) 

where 



(e+e"— Z A ->fc(Z=l))/ 

?('=l) f<A _ ^tot W 

ff (e+e-^^V-)( S ) ' 



Both II(<7 ) and 115(g ) can t> e calculated at large \q | from the same 2-point quark-loop diagrams 
shown in Fig. [M] In the chiral limit, the only difference is in two extra 75 matrices present in IEj,. 
Hence, in perturbative QCD 

U 5 (q 2 )=n(q 2 ), (125) 

at any order in a s . This equation is trivial for the leading-order loop diagrams (Fig. \T%k). Using Dirac 
algebra, it is easy to check that the extra 75 matrices cancel each other (75 = 1) in the absence of masses 
in the propagators. Furthermore, each gluon line inserted in the loop brings two more 7-matrices which 
do not influence that cancellation. 




From d!25t follows the equation of two dispersion relations (I89t and (11221) . To keep the dispersion 
integrals convergent, we differentiate them n times at a certain q 2 < 0. The result is: 



ds ^j = I ds , R{S l ■ (126) 
{s-q 2 ) n J (a-«2V» 



Note that even in the presence of quark masses the corrections to (11251) and to (I126t are very small, 
0(m 2 u d ). If q 2 is not very large, \q 2 \ ~ 1 GeV 2 , the integrals in (11261) are dominated by the contributions 
of the low-mass hadronic states to the corresponding e + e _ — > /i cross sections. In the case of vector 
current, the states are p-meson and its radial excitations. Hence, R(s) represents a resonance curve 
with the peaks located at s = m 2 , m 2 , etc. (see Fig. [T9l . The validity of (11261 ) at arbitrary (but large) 
q 2 < implies that also i?s(s) ~ R(s). We then expect the spectrum of resonances generated by 
the axial- vector current to resemble, in gross features, the p-spectrum. However, experimental data [4| 
reveal a completely different picture. The lowest resonances in the axial-vector channel are the pion with 
J p = 0~ , the axial meson ai(1260) and the radial excitation 7r(1300). Summarizing, there is a clear 
evidence, based on the observed properties of hadrons, that chiral symmetry in QCD is violated much 
stronger than it is expected from Lqcd- 

An additional source of the chiral symmetry violation is provided by the quark condensate. De- 
composing the quark and antiquark fields in (I105t in the left-handed and right-handed components, 

{01^*10) = {OKtiqR + MWiR + MlO) = {0\^ q B^qL+M q R))\0)^0, (127) 

we realize that condensate causes vacuum transitions between quarks of different chiralities. Hence, in 
QCD one encounters a spontaneous broken chiral symmetry, a specific situation when the interaction (in 
this case Lqcd) obeys the symmetry (up to the small 0(m U) d) corrections), whereas the lowest-energy 
state (QCD vacuum) violates it 10 . We conclude that in order to correctly reproduce the properties of 
hadrons and hadronic amplitudes (e.g., correlation functions), one has to take into account the vacuum 
fields, in particular, the quark condensate. 

10 That is quite similar to the electroweak sector of the Standard Model |1|, where the electroweak gauge symmetry is 
spontaneously broken by the nonvanishing vacuum average of the Higgs field. 




(f) 



Fig. 20: Diagrams used to calculate condensate contributions to the correlation functions of two quark currents. Lines with 
crosses denote vacuum fields. 



4.3 Condensate contributions to correlation functions 

Let us return to the correlation functions (II 14l> and dll7t . restoring nonzero u and d quark masses, 
that is, working with full QCD. The vacuum quark-gluon fields generate new contributions to II or 
II5. In addition to the perturbative loop diagrams in Fig. [2] there are diagrams shown in Fig. l20l 
where gluons, quarks and antiquarks penetrate to long distances, being absorbed and emitted by vacuum 
fluctuations. The vacuum fields have characteristic momenta of O(Aqcd)- Therefore, if the momentum 
scale in the correlator is large, Q = \J — q 2 S> Aqcd, one can approximate the vacuum state by a 
set of constant fields. In other words, the virtual quark propagating at short distances/times between 
the points x and 0, cannot "resolve" the long-distance fluctuations of vacuum fields and percepts them 
in a form of averaged static fields. This approximation makes the calculation of diagrams in Fig. |20| 
straightforward. In addition to Feynman rules of perturbative QCD for virtual quarks and gluons, one 
has to form all possible combinations (11031) of vacuum fields and replace them by the corresponding 
condensate densities, so that the 4-momenta of the crossed lines on these diagrams are neglected. For 
example, the product of quark and antiquark vacuum fields in Fig. |20]d has to be replaced by (qq) . More 
details on these calculations can be found, e.g., in the review fi"5l . The result for the vector-current 
correlator has the following schematic form: 

I% 2 ) = IP-V) + E C d (q 2 )(0\O d \0) , (128) 
<f=3,4,... 

where the first term on r.h.s. corresponds to the perturbative diagrams in Fig. El whereas the sum 
contains the contributions of vacuum condensates with dimensions d obtained from the diagrams in 
Fig.|2Ql To compensate the growing dimension of the operators O d , the coefficients C d contain increasing 
powers of 1/Q. Thus, the condensate contributions die away at Q —> 00 and do not alter the perturbative 
asymptotics of the correlator given by W eTt (q 2 ). Furthermore, as soon as we work at relatively large 
Q Aqcd, it is possible to retain only a few first terms in the sum , that is, neglect diagrams with more 
than 3-4 vacuum fields emitted from the virtual quarks in the correlation function. 

For the axial- vector correlation function II5 ((? 2 ) one obtains an expression similar to dl28t with 
the same perturbative part (up to very small corrections of 0{rn\ d )), but with different coefficients C d at 
certain condensate terms. The most important deviation from the vector-current case is in the value and 
sign of C*6, i.e., in the 4-quark condensate terms. Thus, addition of condensate effects leads to an explicit 
violation of (U25t . It is then not surprising that hadron resonances contributing to R(s) and Rs(s) in 
(I126t are different. 



Furthermore, the method of correlation functions allows to reproduce an important relation for 
the pion mass, explaining the smallness of m^. Note that from the point of view of naive quark model 
7r and p mesons differ only by orientations of quark spins. Why is then <C m p and, moreover, 
m-K < Aqcd? We consider the correlation function similar to d 1 1 7I> . but for simplicity, containing 
charged axial currents: 

n%=ij d 4 x e ^(0|T{j M5 (x)ii 5 (0)}|0) = -g^%(q 2 ) + q^s{q 2 ) , (129) 

where = u'-f^^d. This current is the part of the Standard Model weak current d67b . responsible for 
u — > d transition , e.g. the 7r — > fiu^ decay. The hadronic matrix element which determines this decay, 

(0\u W d\ir + (q)) =iq»U, (130) 

is parameterized via the pion decay constant f % which plays an essential role in our analysis. It is 
convenient to multiply ( 11291 ) by q^q^/q 2 , forming a combination of invariant amplitudes: 

""""-nL, = nUg 2 )-5 2 n 5 (g 2 ) = n 5 (g 2 ). 



q 2 

Note that 1I5 = at m U) d = and, in particular, the perturbative part of 1I5 vanishes as 0(m^ d ). Con- 
cerning nonperturbative part, the only first-order in m u ^ contribution is given by the quark condensate 
diagram in Fig.l20h 

ft 2 ) = _ (m u + m d )((uu) + (dd)) + Q(m 2 } (131) 
q z 

To proceed, we use for IT5 the dispersion relation of the type (1122b . To obtain the imaginary part, one 
has to return to the unitarity relation (l75t. identifying |i) with |0) and Ta with the correlation function. 
The result is: 

2Imn 5 ( S ) = E(°l^5|M(^|ii 5 (9)|0) • (132) 

Importantly, only pseudoscalar states contribute to the above sum, because the matrix elements for the 
axial-vector mesons a\ and its excitations vanish, being proportional to the transverse polarization vec- 
tors of these mesons: 

^(OIj^K) ~q%' =0. 
Using (I132t and the definition dl30t we obtain the following expression for the dispersion relation: 

, {-fis6(a-m*)} + ..., (133) 
2 — to > 

where the one-pion state contribution is shown explicitly (with p 2 = m 2 ) including the phase space 
proportional to (^-function. Ellipses denote excited pions and multiparticle states with the same quantum 
numbers. Integrating out <5-function and substituting d 1 3 1 1 > in l.h.s. we obtain 

- (^ + m d )((nu) + (dd)) + 2 = _J^4 _ ^ J>4_ (134) 
q A ' m z — q z '—f mr, — q z 

where the sum over higher states is also shown schematically. To fulfil this equation at large q 2 one has 
to demand that m 2 ~ m u + m^. Simultaneously, the decay constants of excited states have to obey: 



00 

ds 



fa ~ m u + rrtrf, otherwise the q 2 asymptotics of l.h.s and r.h.s. in ( 11341 is different at 0(m u ^) . We 
then reproduce the well known Gell-Mann-Oakes-Renner relation 

- (m u + m d )((uu) + (dd)) + 0{m 2 u4 ) = f 2 m 2 n . (135) 

This relation reflects the special nature of pion in QCD. The anomalously small pion mass is not acci- 
dental and is closely related to the spontaneous chiral symmetry breaking via condensate. If a symmetry 
in quantum field theory is broken spontaneously, there should be massless states (Nambu-Goldstone par- 
ticles), one per each degree of freedom of broken symmetry. The three pions, ir + , tt~ , 7r° play a role 
of massless Nambu-Goldstone particles in QCD. In other words, due to the specific structure of QCD 
vacuum fields, the amount of energy needed to produce a pion state tends to zero. The fact that pions 
still have small nonvanishing masses is due to the explicit violation of chiral symmetry via u, d quark 
masses. 

How large is the quark condensate density? Using m n ± ~ 140 MeV, f n = 131 MeV (4| and 
taking the u, d quark mass values from d3Tl . renormalizing them at Q = 1 GeV, one obtains from (1135b . 
typically, 

{qq)(p = lGeV) ^ (-240 ± 10 MeV) 3 , (136) 

where we assume isospin symmetry for the condensates (qq) ~ (uu) ~ (dd). Not surprisingly, the 
estimated value is in the ballpark of Aqc d ! Being not a measurable physical quantity, the condensate 
density is a scale-dependent parameter. Since r.h.s. of (I135t is determined by hadronic parameters 
/tt, run, which are both scale-independent, the running of the quark condensate should compensate the 
running of the quark mass given in d37l . that is: 

(m)(Q) = (QQKQo)(^ L A ■ (137) 



4.4 Gluon condensate 

The gluon condensate density is another important characteristics of nonperturbative QCD. This param- 
eter cannot be easily estimated from correlation functions with light quarks, because the latter are dom- 
inated by quark condensates. A very useful object, sensitive to the gluon condensate is the correlation 
function of c-quark currents: 

= ij d 4 xe^(0\T{j^(x)f u (0)}\0) = (- M 2 + q^W{q 2 ) , (138) 

where j c = c^^c is the c-quark part of the quark e.m. current j^ m . Following the same derivation as 
in Lecture 3, we write down the dispersion relation for H c (q 2 ) relating Imil c (s) with the ratio R c (s) 
defined as: 

g(e+e- -> charm) 
R c = , t _ T _^ , (139) 

where the cross section a(e + e~ — ► charm) includes hadronic states with cc content produced in e + e~: 
charmonium resonances ( J/ip, ip' , ...), pairs of charmed hadrons etc. In addition to the perturbative 
c-quark loop diagrams in Fig El one has to include also the contribution of the diagrams shown in 
Fig|20b.b.c. with c quarks emitting vacuum gluons. (Remember that nonperturbative quark condensate 
for the heavy c quark is absent.) The resulting relation has the following form: 

oo 

/ 7^^( s ) = n per V,m;?,a s ) + ^(GG)f c (q 2 ,m 2 c ) , (140) 



where is the mass of the lowest J/ip state in this channel. The function UP ert evaluated from the 
massive c -quark loop diagrams has a more complicated form than in the massless case. The function f c 



is the calculable short-distance part of the diagrams in Fig. I20h.b.c. To achieve a better convergence at 
s — > oo, the dispersion relation is usually differentiated n times. Note that in this case the point q 2 = 
is also accessible: c quarks are still highly virtual at q 2 = because the long-distance region starts at 
q 2 ~ Am 2 . The set of power moments obtained from (11401) with 0(a 2 ) accuracy has a form 



ds 



R c (s) 



3C? 2 . 



(Am 



1 + a s (m c )a n + (a s (m c )) 2 a' n + b. 



(ifGG) 



(Am 2 ) 2 



(141) 



with calculable coefficients r n ,a n ,a' n ,b n . The natural scale for a s in the perturbative loops is in this 
case the virtuality Q ~ m c . The moments (I14U are used to extract the gluon condensate density and, 
simultaneously the c quark mass, employing the experimental data on R c (s) on l.h.s. The estimate of the 
condensate density obtained first in [ 17 1 is: 



— (GG) = (330 MeV) 4 ± 50%, 



7T 



(142) 



again within the range of Aqcd- The value of the gluon condensate density is usually given multiplied 
by a s for convenience, because this product is scale-independent. 



5. RELATING QUARKS AND HADRONS: QCD SUM RULES 
5.1 Introducing the method 

The relation (I14U obtained first in [ 16] is a well known example of a QCD sum rule. The method devel- 
oped by Shifman, Vainshtein and Zakharov fTD employs quark-current correlation functions calculating 
them in the spacelike region, including perturbative and condensate contributions. Consider for example 
the correlation function (U14t . with the result of QCD calculation having the form (11281 . Note that in 
the latter expression short- and long-distance effects are separated. The perturbative part W ert (q 2 ) and 
the coefficients Cd(q 2 ) take into account short-distance quark-gluon interactions with characteristic mo- 
menta larger than the scale Q. Both n peri and Cd are process-dependent, i.e., depend on the choice of 
the currents. On the other hand, the condensate densities absorb, in an averaged way, the long-distance 
interactions with momenta less than Q and are process-independent. The universality of condensates 
allows to calculate correlation functions in different channels without introducing new inputs, in an al- 
most model-independent way 1 1 . To obtain the sum rule, the QCD result for the correlation function is 
matched, via dispersion relation, to the sum over hadronic contributions (the integral over hadronic cross 
section). 

A twofold use of the sum rule relations is possible. Firstly, using experimental data one saturates 
the hadronic sum in the dispersion integral and extracts the universal QCD parameters: quark masses, 
a s , condensate densities. One example is the cc sum rule d 14 1 1 > discussed above. Secondly, having fixed 
QCD parameters, one calculates, with a certain accuracy, the masses and decay amplitudes of the lowest 
hadrons entering the dispersion integral. The relation dl35t obtained from the sum rule (I134t can serve 
as an example. Note that dl34t is unique, because all states in the hadronic sum except the pion and all 
QCD terms except the quark condensate are absent at 0(m U: d)- 



5.2 The p-meson decay constant 

To demonstrate in more detail how the method works 12 , let me outline the calculation of the p-meson 
decay constant from the correlation function dl 14I) . The dispersion relation for LT Ml/ can be written in the 

"The expansions similar to i 1281 represent another example of OPE in QCD (see Lecture 3.3), when a product of two 
currents is expanded in a set of local operators Od- The perturbative part in this case is interpreted as a unit operator with zero 
dimension. 

"reviews can be found e.g., in fl5l.fl8l. 



following form: 



^(q) = ^ J g _ g t_ iS \ E(0|iMlM(ftn|>(g)|0) 

a^uQ 2 + q^Qu) \ , (143) 



ds 



s — q 2 — i5 



f -ls(s-ml) + p h (s) * • ■■- 



where the p-meson contribution (h n = p ) is isolated from the sum and the hadronic matrix element is 
substituted: 

<0tf>°> = -^m p e<?>, ( 144 ) 

determined by the p decay constant f p . The integrand p h (s) (spectral density) includes the sum over 
excited and multihadron states. We take into account the experimental fact that the p resonance strongly 
dominates in the low-energy region 2m^ < yfs < 1 GeV, so that p h (s) practically starts from some 
threshold value so ~ 1 GeV. Switching to the invariant amplitude H(q 2 ) and using the result d!28l > of 
QCD calculation, one obtains from (11431) the desired relation, a prototype of QCD sum rule: 



+ J^S - n-(, 2 ) + cV)(o| 0j |o, , ( U5) 

which is valid at sufficiently large \q 2 \. To proceed, one applies to both sides of this equation the Borel 
transformation defined as: 

B M 2U(q 2 )= lim UL± (—s) U(q 2 )^U(M 2 ). (146) 

—q 2 /n=M 2 



This transformation deserves a clarifying comment. Differentiating H(q 2 ) many times in q 2 means, one 
is effectively approaching the long-distance region. Indeed, with an infinite amount of derivatives the 
function H(q 2 ) is defined at any q 2 , including q 2 > 0. The q 2 — > — oo limit works in the opposite 
direction: one penetrates into the deep spacelike asymptotics. Combining the two transformations in 
dl46t at M 2 = —q 2 /n, one fixes the virtuality scale at 0(M). Using the school-textbook definition of 
the exponent: e x = lim„_^ 00 (l + x/n) n , it is an easy exercise to prove that 

B M A^—\=e-<l M \ (147) 

As a result, Borel transformation exponentially suppresses the integral over p h (s) in dl45t with respect 
to the p-meson term. Furthermore, after applying B to r.h.s. of ( 11451 ) the coefficients C<j(M) contain 
powers of l/M 2 . Hence, at large M 2 it is possible to retain only a few low-dimension condensates in 
the sum, e.g., at M 2 ~ 1 GeV 2 a reasonable approximation is to neglect all operators with d > 6. The 
explicit form of the QCD sum rule (I145t after Borel transformation is (T7\ : 

fi e -m*/Ml + J dg p h {s)e s/M 2 
so 



M 2 
2 



1 (-, a a (M) \ | {m u +ra d ){qq) 1 i^G^G^) 112n a a (gq) 2 
i^l 1 " 1 " tt J" 1 W- h 12 W- 81 M b 



(148) 



In obtaining the above relation, the four-quark vacuum densities are factorized into a product of quark 
condensates. The quark-gluon condensate has very small coefficient and is neglected. The running 
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Fig. 21: The p meson decay constant calculated from the sum rule \148\ neglecting all excited and continuum states (solid), 
as a function of the Borel parameter, in comparison with the experimental value (boxes). The dashed curve corresponds to 
an improved calculation , where the sum over excited and continuum states is estimated using quark-hadron duality with a 
thresholds 1 ^, = 1.7GeV 2 . 




coupling a s is taken at the scale M i.e., at the characteristic virtuality of the loop diagrams after Borel 
transformation. 

Importantly, there exists a SVZ region [ 17 1 of intermediate M 2 where the p meson contribution 
alone saturates the l.h.s. of the sum rule A148I ). To illustrate this statement numerically, in Fig. l2~llthe 
experimentally measured f p (obtained from the p° — > e + e~ width) is compared with the same hadronic 
parameter calculated from the sum rule (I148t where all contributions of excited and continuum states are 
neglected. One indeed observes a good agreement in the region M 2 ~ 1 GeV 2 . 

An important step to improve the sum rule (11481) is to use the quark-hadron duality approximation. 
The perturbative contribution to the correlation function (the sum of Fig.EJdiagrams) is represented in 
the form of a dispersion integral splitted into two parts: 

7*/^(.) fft^M (149) 



s — q 2 — iS J s — q 2 — id 

so 

The integral over the spectral function p (s) in (11451) is approximated by the second integral over the 
perturbative spectral density p per *(s) in dl49t . The latter integral is then subtracted from both parts of 
(U45t . Correspondingly (I148t is modified: the l.h.s. contains only the p term, and, on the r.h.s., the 
perturbative contribution has to be multiplied by a factor (1 — e~ s °' M ). The numerical result obtained 
from the duality-improved sum rule (I148t is also shown in Fig.|^ The agreement between the sum rule 
prediction and experiment is impressive: 

j(QCDSR) = 213 MeV ± ( 10 _ 15)% (150) 

whereas f p xp = 216 db 5 MeV |4|. The estimated theoretical uncertainty quoted in (11501) is typical for 
QCD sum rules, reflecting the approximate nature of this method. So far we only reproduced f p , using 
the experimental value of m p . In principle, it is possible to go one step further and estimate also the 
p-mass from the sum rule . Furthermore, a very similar sum-rule analysis for the correlation function of 
axial currents ( 11291 ) (the invariant amplitude IT5) successfully reproduces the value of f n (see Fig.l2^t. 




Fig. 22: The tt meson decay constant calculated from QCD sum rule [171 for the correlation function of axial-vector currents ( 
solid line), in comparison with the experimental value (boxes). The quark-hadron duality threshold is fitted simultaneously as 
Sq — 0.7 GeV 2 . The uncertainty of about 10-15% has to be added to the theoretical curve. 



5.3 Baryons 

It is possible to extend the method of QCD sum rules to baryons |[T9l l20l . The idea is to construct 
special quark currents with baryon quantum numbers which can serve as a source of baryon produc- 
tion/annihilation from/to QCD vacuum. In reality, such currents hardly exist 13 , but they are allowed in 
QCD if colour-neutrality is obeyed. A well known example is the Ioffe current with the nucleon quantum 
numbers (i.e., spin 1/2): 

J N (x) = e abc (u aT (x)C^u\x)h 5 ^d c (x) , (151) 

where a, b, c are color indices, C = 7270 is the charge conjugation matrix and, for definiteness, the 
proton flavour content is chosen. As a next step, one constructs a correlation function 

U N {q) = i J d 4 xe iqx (0\T{J N (x)J N {0)}\0) . (152) 

The corresponding diagrams are shown in Fig.|23j including perturbative loops and vacuum condensates. 
They look quite different from the quark-antiquark loops, but differences concern the short-distance parts 
of the diagrams. The universality of condensates allows one to calculate Hn(q) without introducing new 
input parameters. The hadronic contribution contains a total sum over states produced and annihilated 
by the current J , starting from the lowest possible state, the nucleon: 

(0\J N \N)(0\J N \N) 

H-n(q) = 9 5 1- {excited resonances, multiparticle states} . (153) 

mfj - q l 

The derivation of the QCD sum rule is done along the same lines as in the previous subsection. Omitting 
the details, let me only mention that from this sum rule an approximate formula for the nucleon mass is 
obtained, 

m N ~ [— (2.0)(2tt) 2 <0|^|0)(^ = lGeV)] 1 / 3 , (154) 

13 In models of Grand Unification predicting proton decay via intermediate superheavy particles, the currents we are dis- 
cussing are realized effectively in a form of localized 3-quark operators annihilating the proton. 
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Fig. 23: Some of the diagrams contributing to the correlation function of two baryon currents: (a) the lowest-order tree-quark 
two-loop diagram; (b) the 0(a a ) correction, (c),(d) -quark condensate diagrams. 

relating it to the quark-condensate density. In fact, the quark masses m u ^ themselves generate very 
small corrections to this relation and are neglected. Thus, QCD sum rules provide an answer to the 
question that was raised in Lecture 1: almost 99% of the baryonic mass in the Universe are due to 
vacuum condensates. 

5.4 Quark mass determination 

As already mentioned, the sum rule (11411) can be used to extract the c-quark mass. Here, the role of Borel 
transformation is played by a simple differentiation at q 2 = 0, which turns out to be more useful. At 
low n, the moments d!41t are especially convenient for m c determination because the gluon condensate 
effects are small. Replacing c — > b, tp — > T , etc., one obtains analogous sum rule relations for b quark, 
where gluon condensate is much less important, being suppressed by m^ 4 . Recent analysis [21 ] of these 
sum rules yields m c (m c ) = 1.304 ± 0.027 GeV, m b (m b ) = 4.209 ± 0.05 GeV. Another subset of 
charmonium sum rules (higher moments at fixed large q 2 < 0) was recently employed in [22|, with a 
prediction for m c in agreement with the above. 

The heavy-quark mass determination using sum rules is also done in a different way, employing 
the large n moments (I14H which are less sensitive to the cross-section above the open flavour threshold. 
These moments, however, demand a careful treatment of Coulomb interactions between heavy quark 
and antiquark in the perturbative diagrams. Remember that one-gluon exchange yields Coulomb po- 
tential. Close to the threshold of heavy quark-pair production, yfs ~ 2rriQ ( Q = b,c) this part of 
quark-antiquark interaction becomes important, and at large n the near-threshold region dominates in the 
perturbative coefficients a n , a' n in (I14U . A systematic treatment of this problem is possible in nonrela- 
tivistic QCD (NRQCD), a specially designed effective theory obtained from QCD in the infinite heavy 
quark mass limit (for a review see, e.g. [23 1). 

Equally well, QCD sum rules allow to estimate the masses of light u, d, s quarks. To give only one 
typical example from the vast literature let me refer to the recent analysis l24l . based on the correlation 
function for derivatives of the s-flavoured vector current j SfJL = sj^q, (q = u, d): 




(155) 



The QCD answer for II s is proportional to (m s — m q ) 2 ~ m 2 (due to d^s = ( m s — r nq)sq), making this 
correlator very sensitive to m s . One calculates the usual set of diagrams shown in Figs.[W]and|20l where 
the quark lines are now s and q. Furthermore, the recent progress in the multiloop QCD calculations 
allows to reach the O(a^) accuracy in the perturbative part of II s and to include also 0(a s ) corrections 
to the condensate contributions. The hadronic spectral density pk{s) for this correlator is saturated by 
J p = + states with s -flavour, e.g. Ktt states with L = 0. The sum rule has the form 

J dsp K (s)e~ s/M2 = (m s - m u ^) 2 

Data of kaon S'-wave scattering on n, rj, rj were used to reproduce p K (s). The resulting prediction [24 1 
for the mass is 

m s (2 GeV) = 99 ± 16 MeV. (157) 

The ratios of the light (u, d, s) quark masses can be predicted in QCD from the relations for pions 
and kaons, similar to (U35t . A systematic derivation is done employing chiral perturbation theory, an 
effective theory obtained from QCD in the a low-energy limit, using instead of quarks and gluons, the 
pion and kaon degrees of freedom (for a review see [25 1). The result is [26 1: 

Tfl Tfi ^TTl 

— = 0.553 ± 0.043, — = 18.9 ±0.8, — = 24.4 ±1.5. (158) 

m d m d m u + m d 

From the above ratios and the value d!57t one obtains m u (2 GeV) = 2.9 ± 0.6 MeV and m d (2 GeV) = 
5.2 ± 0.9 MeV. 

We see that QCD sum rules are extremely useful for the quark mass determination. The m q 
values extracted from sum rules are included, together with the lattice determinations and results of other 
methods, in the world-average intervals in [4| presented in d3TT >. 



UP ert (M) + + 



(156) 



5.5 Calculation of the B meson decay constant 

In 5-meson decays the CKM parameters of the Standard Model are inseparable from hadronic matrix 
elements. Hence, without QCD calculation of these matrix elements with an estimated accuracy, it is 
impossible to use experimental data on B decays for extracting the Standard Model parameters and for 
detecting/constraining new physics effects. Currently, lattice QCD provides many hadronic parameters 
for B physics, with a continuously improving accuracy. QCD sum rules represent another actively used 
working tool. With condensates and quark masses determined from a set of experimentally proven sum 
rules for light-quark and heavy quarkonium systems, one has a real possibility to assess the theoretical 
accuracy of the sum rule predictions by varying the input within allowed intervals. 

One of the most important applications of QCD sum rules is the determination of the S-meson 
decay constant fg defined via the matrix element 

m b (0 | qi lh b | B) = m 2 B f B , (159) 

(q = u, d, s). Note that fg multiplied by V U b determines the width of leptonic B decays, such as 
B~ — ► t~v t . To calculate fs from QCD sum rules, one usually employs the correlation function: 

nf \q) = i J d 4 xe^^0\T{m b q(x)i l5 b(x),m b b(0)i 75 q(0)}\0) , (160) 

so that the lowest B meson term in the hadronic sum for this correlation function contains 



To obtain the sum rule, one needs to calculate H5 (q) from the perturbative and condensate diagrams 
in Figs. [2] and with 75 vertices emitting and absorbing b and q lines. One, therefore needs an 
input value for m& and m s (m Uj d can safely be neglected). These values are taken from the analyses 
overviewed in the previous subsection. The sum rule has a form similar to d!48t for f p , Naturally, 
the expressions for the perturbative part and for the coefficients are completely different. Also the 
hierarchy of contributions in the heavy-light correlation function differs from the light-quark case. Now 
the quark condensate becomes very important being proportional to m^qq). The recent updates of the 
sum rule for Jb obtained in l27l l28l take into account the 0{a 2 s ) corrections to the heavy-light loop 
calculated in [29 1. The numerical prediction of the sum rule, taking m^m?,) = 4.21 ± 0.05 GeV, is: 
fs = 210 ± 19 MeV and f Bs = 244 ± 21 MeV j27]], in a good agreement with the most recent lattice 
QCD determinations. 

5.6 Light-cone sum rules and B —* tt form factor 

To complete our brief survey of QCD sum rules, let me introduce one important version of this method, 
the light-cone sum rules (LCSR) 11301 071 used to calculate various hadronic amplitudes relevant for 
exclusive processes. In the following, we consider the application of LCSR to the B —> it transition 
amplitude (see Fig. |9j>. The latter is determined by the hadronic matrix element 



generated by the b — ► u weak current <67b . Due to spin-parity conservation only the vector part of the 
current contributes. There are two independent 4-momenta p and q, and one independent invariant q 2 , 
the momentum transfer squared. The initial and final mesons are on shell, p 2 = m 2 and (p + q) 2 = m 2 B . 
It is quite obvious that one needs two invariant functions of q 2 , the. form factors f^{q 2 ) and fg^iq 2 ), to 
parameterize this matrix element. Only one form factor f£ is interesting, the other one is kinematically 
suppressed in the measurable B — ► nlvi semileptonic decay rate (I = e, /x). 

To derive LCSR for f Bn (q 2 ) one uses a specific correlation function, which itself represents a 
hadronic matrix element. It is constructed from the product of the weak uj^b current and the current 
rribbi^d used to generate B in dl59t . The currents are taken at two different 4-points and sandwiched 
between vacuum and the one-pion state. The formal definition of this correlation function reads: 



where it is sufficient to consider only one invariant amplitude F. Note that we now have two kinematical 
invariants, q 2 and (p + q) 2 , still p 2 = m 2 . 

Diagrammatically, the correlation function dl63t is represented in Fig. [24] At spacelike (p + q) 2 < 
and q 2 < 0, very similar diagrams describe the process 7*7* — > ir° , the one-pion production by two 
virtual photons via e.m. currents 14 . One only has to replace all quarks in the diagrams in Fig. |24] 
by either u or d quarks. Both objects: the heavy-light correlation function (I163t and the 7*7* — > it 
amplitude contain one virtual quark propagating between vertices and a quark-antiquark pair which is 
emitted at points x and and converted into a real pion state. At large spacelike external momenta 
\{p + q) 2 \, \q 2 \ 3> Aq C£) the space-time interval x 2 ~ approaches the light-cone. Hence, the virtual 
quark in both (two-currents — ► pion) amplitudes propagates at short distances allowing a perturbative 
QCD description. The calculable short-distance parts are process-dependent. In case of the correlation 
function (I163t we have a virtual b quark propagating between vertices of flavour-changing currents, 
whereas in the 7*7* — ► 7r° amplitude the light quark propagates between e.m. vertices. The long- 
distance part in both cases is, however, the same vacuum-to-pion matrix element of light quark and 

14 This process is experimentally accessible in e + e~ —* e + e~ir° two-photon (double-tagged) collisions. 



< 7T + (p)\u ltl b\B(p + q) >= 2f+ n (q 2 ) Pfl + f^(q 2 ) + f^(q 2 ) q, 



(162) 




F(q 2 , (p + qf)p^ + F(q 2 , (p + q) 2 )q p , 



(163) 



p + q 



<J 



Fig. 24: Diagrams contributing to the correlation function \163\ . from left to right: the leading-order, 0{a s ) correction, soft 
gluon. The blob with ty denotes the pion distribution amplitude. 



antiquark emitted at points x and 0: 

<7r(p) I q X {x)T a q 2 {0) | 0) , (164) 

where in one case q\ = u, q 2 = d, it = tt + and in the other case q\ = u(d), q 2 = u(d), it = ir°. Due to 
isospin symmetry, the difference between these two configurations is indeed very small. In order to treat 
the diagram with one extra gluon entering the pion, one has to introduce an additional quark-antiquark- 
gluon matrix elements of the type 

<7r(p) I u(x)g s G^(y)T b d(0) | 0) , (165) 

where x 2 ~ y 2 ~ (x — y) 2 — > 0. In the above T a j, denote certain combinations of Dirac matrices. 

Having separated short- and long-distances, one is able to calculate the correlation function (1163b 
in a form of light-cone OPE, where the short-distance part (the virtual b quark propagator plus gluon 
corrections) is multiplied by a long-distance part, the universal matrix elements such as dl64t . dl65t . 
The latter can be parameterized in terms of light-cone distribution amplitudes of the pion [32], The most 
important of them is defined by 

(7T(p)\u(x) lfll5 d(0)\0) = -ippU [ due^^u,^), (166) 

J 

where \i is a characteristic momentum scale, determined by the average x 2 in the correlation function. 
In the above definition, u and 1 — u are the fractions of the pion momentum p carried by the quark 
and antiquark, in the approximation where one neglects the transverse momenta of the constituents with 
respect to the longitudinal momenta. The leading-order answer for F obtained from the first diagram in 
Fig-GHusing (|166t is quite simple: 



du (p n (u, fl) 
m 2 — (q + upY z 



F{q 2 , (p + q) 2 ) = m b U I 7 ,1,, ■ d67) 



The pion DA ip n (u) plays here the same role of nonperturbative input as the condensates in the conven- 
tional QCD sum rules considered above. Asymptotically, that is at fi —> oo, QCD perturbation theory 
implies 99^(^,00) = 6u(l — u). However, at the physical scale fi ~ m&, at which the OPE is ap- 
plied to the correlation function (11631) . nonasymptotic effects also contribute, which we will not discuss 
for brevity. Importantly, there is a power hierarchy of different contributions stemming from the dia- 
grams in Fig. I231 determined by the large scale in the correlation function. This scale is given by the 



p + q 



p + q 



Fig. 25: Diagramatic representation of the dispersion relation for the correlation function M63\ ) 



virtuality of b quark, i.e. by the quantity which stands in the denominator of the b -quark propagator: 
m b ~ (l + u p) 2 = m b ~ (P + l) 2u ~ 1 2 {^ ~ u )- Importantly, this quantity remains large when q 2 
is positive (timelike) but not very large, q 2 -C m 2 , allowing one to penetrate into the lower part of the 
kinematical region < q 2 < (tub — ni^) 2 of B — ► irlui decay. An example of a subleading contribution 
is the diagram with an additional gluon entering pion DA (the third one in Fig. l24t . It has two powers 
of inverse scale and is suppressed with respect to the leading-order diagram. Over recent years 0(a s ) 
corrections and quark-antiquark-gluon contributions have been calculated to improve the result (I167t . 

Having evaluated F as a function of q 2 and (p + q) 2 , we are still half-way from the final sum 
rule. The next important step is writing down a dispersion relation in the variable (p + q) 2 , the external 
momentum of the current with B meson quantum numbers. The diagrammatical representation of the 
dispersion relation is shown in Fig|25] It contains the same set of hadronic states as in the sum rule for 
fs, starting from B meson ground state. Using (I159E and (I162t one obtains the hadronic representation 
we need: 

,2 (m , „^2^ _ 2m2 BfBft^P 2 

mb{m 2 B — (p + q)' 

The ellipses in the above denote the contributions from the excited B and from continuum states. Equat- 
ing the QCD result (11671) in the region of validity (p + q) 2 < with the dispersion relation d 168l > one 
obtains a raw sum rule relation for fsfs (p 2 )- The rest of the calculation follows the usual QCD sum rule 
procedure: Borel transformation in (p + q) 2 and subtraction of the contribution from higher states invok- 
ing quark-hadron duality. One finally arrives at an expression for the desired form factor, the dominant 
term of which obtained directly from dl67t is given by 

<•+ / ?\ firm 2 f 1 du , . . . / m 2 n m? — q 2 (1 — u)\ m 
ft(<l 2 ) = ^- J k- / — M",M6) + .--)e*P -§- b XA= '- . (169) 



F(q 2 , (p + q?) = J',*?:,'™ +■■■■ ( 16§ ) 



Here, M is the Borel mass parameter, and the scale //& reflects the characteristic virtuality of the correla- 
tion function, n 2 = m 2 B — m\. The integration limit A = (m\ — p 2 ) / '(sq — p 2 ) depends on the effective 
threshold Sq above which the contribution from higher states to the dispersion relation dl68t is canceled 
against the corresponding piece in the QCD representation (1167b . The parameters fs,s B are usually 
taken from the sum rule for considered in the previous subsection. The most recent predictions for 
IbM 2 ) obtained from LCSR [33] were used to extract | V u b \ from the measurements of the exclusive 



semileptonic width at B factories l34ll . using the formula for the decay rate: 



dT(B -> ttW) _ G 2 \V ub 
dip ~ 24vr 3 



(El - m: 



l) i/2 [ftM 2 )\ 2 ■ 



(170) 



where E n is the pion energy in the B meson rest frame. 
6. CONCLUSIONS 

QCD has a thirty-years history and embraces many approaches, some of them developed quite indepen- 
dently from the others. Due to self-interactions of gluons this theory has an extremely rich dynamics, 
combining asymptotic freedom at short distances with the self-emerging energy scale Aqcd and con- 
finement at long distances. Accordingly, QCD has two different phases: the perturbative one responsible 
for quark-gluon processes at large momentum transfers and the nonperturbative one where the only ob- 
servable states are hadrons formed by confined quarks and gluons. Yet there is no complete analytical 
solution for the hadronic phase of QCD. One has to rely on approximations: either numerical (QCD 
on the lattice) or analytical (QCD sum rules). In addition, several effective theories corresponding to 
different limits of QCD and exploiting the rich symmetry pattern of the theory are successfully used. In 
the overview of QCD given in these lectures I tried to emphasize the importance of the hadronic aspects 
of QCD, where the most nontrivial phenomena and challenging problems are accumulated. 
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